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Abstract. Any subset of the plane can be approximated by a set of square pixels. This transition 
from a shape to its pixelation is rather brutal since it destroys geometric and topological information 
about the shape. Using a technique inspired by Morse Theory, we algorithmically produce a PL approx- 
imation of the original shape using only information from its pixelation. This approximation converges 
to the original shape in a very strong sense: as the size of the pixels goes to zero we can recover impor- 
tant geometric and topological invariants of the original shape such as Betti numbers, area, perimeter 
and curvature measures. 



A common problem in computational topology is to try to recover an object embedded in M" when 
only some distorted version of it is known. Inspired by digital imaging we consider specifically the 
problem of a pixelated subset of the plane, which is to say that it has been replaced by a set of square 
pixels on a grid. Since it is common to represent images by a grid of small pixels, it seems intuitive 
that any object could be recovered from its pixelations. Indeed, when the resolution of an image is 
fine it is difficult to detect any difference from the original using only the human eye, and we may 
expect to be able to recover even deep geometric invariants.lated subset of the plane, which is to say 
that it has been replaced by a set of square pixels on a grid. Since it is common to represent images 
by a grid of small pixels, it seems intuitive that any object could be recovered from its pixelations. 
Indeed, when the resolution of an image is fine it is difficult to detect any difference from the original 
using only the human eye, and we may expect to be able to recover even deep geometric invariants. 

A pixelation associated to a subset of a plane is simply the set of all pixels in a grid of a certain 
size which touch that subset (see Figure [T]l. By changing the size of the grid, we can view pixelations 
of the object with finer and courser resolutions. We call the pixelation associated to a set 5 C 
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created from a grid of side-length e Pe(S'), or tlie e-pixelation of S (see Definition [lT]). The variable 
e is called the resolution of the pixelation. 




Figure 1 . An example pixelation of a line with slope ^ 

For small resolutions these pixelations seem to resemble the original shape. However, we need a 
precise notion of what it means to resemble the original object. At the very least it seems that the 
important topological and geometric invariants of the pixelation will be close to the corresponding 
variants of the original set. Even when viewing simple examples, we can quickly see that not all 
invariants will converge for small resolutions. Most obviously, the total curvature of Pe{S) will tend 
to increase to infinity, since the pixelation contains many comers. Similarly length of the boundary 
may not converge. To see this consider the line segment from (0, 0) to (1,1). The length of the line 
is y/2, but each of its pixelations will have a boundary whose length is close to 4. 

We see that pixelations do destroy various geometric invariants. However, it may be that con- 
vergence of these invariants is too much to ask for. A much simpler request is for the pixelation to 
eventually converge in Euler characteristic. Unfortunately not even this is guaranteed. Some pix- 
elations will contain fake cycles which do not correspond to any cycle of the original set (indeed 
arbitrarily many may appear even if the original set was contractible). Worse yet, these fake cycles 
may not disappear for small resolutions. Section|A]has further information on this phenomenon, but 
for now it is sufficient to consider the example of a line of slope 1 and a line of slope | (which is 
shown in Figure [2]). The e-pixelations have certain scale invariance properties which guarantee that a 
fake cycle will always appear, even for very small e. 

We see that the process of pixelation can destroy a great deal of mathematical information of the 
set. However, to the eye, the pixelation can look very close to the original set. Therefore it seems that 
it should be possible to use only the pixelation to reconstruct the original set. In this paper we will 
create an explicit algorithm which can reconstruct a PL set S from only its pixelations. However, to 
motivate this algorithm and to generate the theorems needed to show that this algorithm works, we 
will need to examine simpler examples of pixelations. 

In Section[2]we deal with pixelations of the graphs of functions. Theorem 2.2 an analogue of 



the intermediate value theorem, tells us that these pixelations must be contractible. Therefore we do 
not have to worry at this stage about fake cycles. To approximate the graph of a function we attempt 
to use a secant line approximation. This is done by connecting points within the pixelation of the 
function to create a PL approximation. However, additionally error will be generated from the fact 
that our sample points are being selected not from the graph of function, but merely near to the graph. 
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Figure 2. The pixelation of the angle A (l, |) contains two holes. 



Theorem 2.6 gives us bounds on the error between a secant Une of our function and the approximate 
secant Une from the pixelation. Note that these bounds depend both on the resolution e and the width 
of the secant line. We want to choose our PL approximation in a strict algorithmic way such that these 
errors vanish. To do this we use the notion of a spread function. 

A spread function is simply a function a : — Z+ such that the following conditions hold: 

lim a{e) = oo, lim ecr(e) = 0. 

The spread function tells how wide the secant lines in our approximations should be in terms of the 
resolution of the pixelation. That is, since a pixel is defined to be e wide, the secant lines should 
be arranged to be about ea{e) wide. Theorem 2.11 shows that the two conditions of the spread a 
guarantee Sobolev W^'^ convergence of the secant line approximation. 

In Section |3] we generalize the approximation technique to handle elementary regions, which are 
defined to be regions which line in between the graphs of two functions (note that the graph of 
a function is just a special case of an elementary region). In this case we also cannot have fake 
cycles. We can easily generalize the secant line approximation of functions to elementary regions by 
approximating the top and bottom regions of the elementary region. It is then simple to show that we 
maintain the previous Sobolev VF^'P convergence on the boundary of the approximation. 

Here the notion of a spread becomes very useful. By restricting the spread function a, we can 
dramatically improve the convergence of the approximation. Indeed if require that 

lim e{a{e))'^ = oo, 

so that a increases very quickly, but still does not increase faster than e decreases, we can recover the 
total curvature of the boundary of the elementary region. This is shown in Propositions 3.4 3.6 and 
Corollary 3.7 for sets of increasing generality. However each of these sets is defined as the region 
between the graphs of two functions. 

At this point although we have created an algorithm which strongly approximates an elementary 
regions, but we have completely ignored simple sets such as the union of two intersecting lines. Such 
a set is potentially difficult to approximate, since its pixelation may have more cycles than the original 
set. In Section|4] we tackle the problem of homotopy type as it relates to the pixelations of PL sets. 

Corollary |4.2| ofTheorem |4~T] shows that eventually there is a bijective correspondence between the 
connected components of the pixelation and the connected components of the original set S. A cycle 
in a (reasonably behaved) planar set corresponds to a hole, i.e., a bounded connected component of 
the complement of the set. We can again use Theorem |4.1| to show that every hole of the original set 
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S will eventually correspond to a hole of the pixelation. Therefore any defect in homotopy type is 
caused by the addition of holes in the process of pixelation. This means that to ensure convergence in 
homotopy type, we need only to delete the extra holes using only information from the pixelation. 

An intuitive way to distinguish fake holes from real holes is to note that any real hole must take up 
an actual area in the plane. However, fake holes tend to be composed of a relatively small number of 
really small pixels. Therefore for small resolutions e the fake cycles will have small areas. Appealing 
as it may sound, this idea is difficult to implement rigorously because we do not have an accurate way 
of defining what "small area" means. We rely instead on a more robust approach inspired from Morse 
theory. 

We consider the linear function on the Cartesian plane that associates to each point its z-coordinate. 
For simplicity we assume that its restriction £3 to our PL set S is a stratified Morse function in the 
sense of Goresky-MacPherson, |l9l. In our case this simply means that no two vertices of our PL set 
S lie on the same vertical line. 

The topology of the level sets of £3 is determined by the counting function n 

n{x) := the number of components of £^^{x). 

A pixelated version of n is the function n^, where ne(x) is the number of connected components of 
the column of Pe{S) located at xq. Note that for a cycle to appear in S, the function ns must vary 
and likewise with the pixelation and n^. Therefore determining whether a cycle in Pe{S) is a fake or 
really corresponds to a cycle in S" is a matter of determining how closely the function agrees with 
ns. 

The Separation Theorem (Theorem 4.4 1 states that if x is not a critical value of £3 then n(x) = 
n£{x) for all e sufficiently small. The critical values of £3 correspond to points jumping points of n, 
i.e., points of of discontinuity of n. Using the Separation Theorem we prove several results indicating 
that the jumping points of and the jumping points of n5 are not far apart. This culminates in 
Proposition 4. 10 which states that in the PL case jumping points of must appear close to jumping 
points of and vice versa. Here "close" means within eu^S) where i^(-S') is an integer determined 
by S and called the noise range. 

We ultimately want to classify an interval around every jumping point of as "noise." To do this 
we must estimate the noise range i^{S), but this is determined by the original set S and thus unknown 
given only the pixelation. However if we estimate the noise range using a spread a note that for small 
e we can guarantee that a{e) > i^{S) while ea{e) is small. Therefore if we classify an interval of 
width about ea{e) around each jumping point of as noise, we can ensure that eventually outside 
of the noise = ns while the noise region remains small. 

The noise intervals are important since they are chosen to eventually contain all fake cycles of 
Pe{S) while taking up a vanishingly small portion of the real line. Therefore to approximate homo- 
topy type within noise intervals we need only cover every cycle (since any cycle appearing will be a 
fake cycle). Since noise does not take up much of the plane, we do not need to be careful as to how 
we approximate within noise intervals, so we will do this by covering each connected component of 
the noise with the smallest rectangle which covers it. 

Call the components of the complement of the noise intervals regular intervals. By the definition 
of noise, the regular intervals contain no jumping points of n^. Therefore is continuous on the 
regular intervals, which implies that the parts of PeiS) which lie over regular intervals look like strips 
that span the entire regular interval. Each strip can be interpreted as the pixelation of an elementary 
set and we have explained how to deal with such objects. 

Putting together all these facts we obtain Algorithm 4.12 which associates to each e-pixelation a 
PL set. This algorithm works by using the function n^ and the spread a to divide the pixelation into 
noise and regular intervals, and then approximates within each interval using the appropriate results. 
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That is to say within noise intervals it covers each connected component with a rectangle and within 
regular intervals it connects the tops and bottoms of every a-th column. An example of the result of 
this algorithm can be seen in Figure [3] and the Algorithm is restated in terms of a computer program 
in Appendix PI 



Figure 3. The result of Algorithm 4.12 as applied to the intersection of two red 



lines. Note that near the point of intersection the approximation resembles a rectan- 
gle, since this lies within a noise interval. 



We wish to state that Algorithm 4. 12 captures all the important geometric and topological invari- 
ants. To do this we use the language of normal cycles. These are 1 -currents with supports contained 
in the bundle of unit tangent vectors. Appendices [B] and [C] are included for readers unfamiliar with 
these objects. For purposes of this general discussion it suffices to mention that the normal cycle of 
a subanalytic subset S of the plane is a piecewise closed curve in the unit tangent sphere bundle 
object which encodes geometric and topological information about of the set S through integration of 
certain canonical 1-forms. 

In general the normal cycle can be thought to be approximated by the collection of all unit outer 
normal vectors of the set. For example, the normal cycle of a bounded domain with C^-boundary 
is the graph of the Gauss map of the boundary. The normal cycle of the square resembles Figure [4] 
Of course the normal cycle is actually a subset of the sphere bundle of the plane, not the plane itself. 
The fiber of the normal cycle over a point in the interior of an edge each point will consist of only 
a single outer normal vector. The fiber over a corner consists of an entire quarter-circle (all possible 
outer normal vectors for the comer). Therefore the length of the normal cycle will be the perimeter 
of the square plus 2-k, or the perimeter of the square plus its total curvature. 

If we prove that the normal cycles of our approximations converge weakly to the normal cycle 
of the original set, we will have shown that our approximation recovers the important geometric 
information of the set. To prove this we use the powerful General Convergence Theorem proved by 
Joseph Fu in 0. This result is restated in this paper as Theorem |5.5| This theorem tells us that 



to have convergence in normal cycles we must prove three things about our approximations. First, 



they must all fit within some compact subset of the plane. This is clear from the Algorithm 4.12 
Secondly, the mass of the normal cycles of the approximations cannot explode. Since the mass of 
the normal cycle is highly dependent on the total curvature and the perimeter of the set, this follows 
from Corollary |3.7| Finally we must have the Euler characteristic of the approximation converge to 
the Euler characteristic of the original set when restricted to an arbitrary half plane. This requirement 
of the approximation theorem is the most challenging, but it turns to be true as well. Therefore Fu's 
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Figure 4. A representation of the normal cycle of the square. 
General Convergence Theorem implies convergence of our approximations in normal cycle (Theorem 



5.11. 



Theorem 5.1 is a powerful result with many corollaries, since the normal cycle encodes many 



important invariants. Most particularly it implies that the approximation generated by Algorithm 



4.12 will recover both the Euler Characteristic as well as the total curvature on the boundary of any 



generic PL set. 
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1. Definition OF PIXELATIONS 



In this section we give the definition of a pixelation and detail the most basic of properties of a 
pixelation. Questions about recovery of the original set are not tackled until future sections. 

Definition 1.1. (a) Let e > 0. Then we define an e-pixel to be the square in of the form 



The e grid of pixels is the collection of all such e-pixels. 

(b) A union of finitely many e-pixels is called an e-pixelation. The variable e is called the resolution 
of the pixelation. 

(c) For any bounded subset 5 C we define the e-pixelation ofS to be the union of all the e-pixels 
that intersect S. We denote the e-pixelation of S by P£{S). The pixelation of a function / is defined 
to be the piexelation of its graph r(/). We will denote this pixelation by P£(/). 

The variable e is called the resolution, because we think of a pixelation as a representation of 
a computerized image. A smaller choice of e will cause a shape to be approximated by a greater 
number of pixels, which is like saving an image as a higher resolution file. 

Observe that if || • ||oo : [0, oo) is the norm 



then the the pixel Sij{e) can be identified with the || • ||oo -closed ball of radius e/2 and center 



Therefore a pixelation can be thought of as set of points chosen near the original shape. However, 
rather than being chosen randomly these points are chosen to be the closest points on a regular lattice. 
This by no means makes the problem of a recovering a set from its pixelation trivial. As we will see, 
even the homotopy type can be lost during the process of pixelation. Still, the added structure of a 
pixelation will allow us to make stronger approximations than can be done with regular pixelations. 
For example it is reasonable to expect geometric information like total curvature of the boundary to 
be preserved in an approximation. 

The e-pixelation associated to a bounded set S is the collection of pixels in this grid that contain 
a part of S. A pixelation is a sort of "fattening" of the set, since we include a pixel if any part 
of the set lies within the pixel (see Figure |5]) In this way it is somewhat similar to a tube around 
the set. However a pixelation differs from a tube in that the boundary of any pixelation will have 
corners, and that the distance from the boundary to the set may vary significantly (since pixelations 
are formed by taking a collection of squares in a grid.) As we will see in the following sections, the 
fact that pixelations approximate via a grid means that the geometric properties of pixelations may 
not converge to the geometric properties of the original set as e — (for example, the normal cycle 
of the pixelation will not converge to the normal cycle of the original set.) The good news is that 
geometric information about the original set can be obtained from its pixelations in indirect ways. 

It is easy to see that the pixelation of an set can be fundamentally different from a set. For example, 
every pixelation of a line segment will have many comers (as long as the line is not parallel to the x or 
y-axis.) This impUes that as e becomes small, the total curvature of the boundaries of pixelations of 
a line segment will increase in an unbounded fashion. This contradicts the fact that the line itself has 
no curvature and shows that essential geometric features of the set are not preserved in the pixelation. 

In fact, though it is more difficult to see, pixelations do not preserve the topological properties of 
the set. Consider the set S consisting of a line of slope ^ and a line of slope 1 which both start at the 



Il(^,y)ll 



max{|x|, |y|} 
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Figure 5 . An example pixelation of a line with slope ^ 



origin. Then it can be shown that for every e, Pe{S) contains a hole which prevents Pe{S) from being 
contractible (see Figure [6]) The position of the hole depends on e, but it will always exist for any e. 
This example is far from unique. In fact, by altering the slope of the two lines we can create a set 
whose pixelations contain any desired number of holes. This shows that pixelations do not preserve 
the topological properties of the underlying set. 




Figure 6. If S consists of a line of slope | and a line of slope 1 starting at the 
origin, then Ps{S) is not contractible for any e. 
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However, when e is very small, it is easy for the human eye to detect the kind of shape that 
probably generated a pixelation. Mathematically this corresponds to the fact that it is possible to 
create an approximation from the information in the pixelation which will respect topological and 
geometric features. The approximation can be algorithmically generated, but we delay defining it 
until much later in the paper (since the motivation for it rests on many properties of pixelations.) 

The remainder of the paper will examine specifically what pixelations preserve and what they 
destroy. Furthermore it will explain how to recover the information "lost" in the process of creating a 
pixelation. 

We end this section with a few basic properties of pixelations which will be useful. 
Proposition 1.2. For any two sets U and V we have 

P,{UUV) = Pe{U)UPe{V) 



Proof. This follows immediately from the definition of a pixelation as the union all pixels which 
intersect a set. □ 

Proposition 1.3. For each e define T^{S) to be 

T^{S) := {x eR^ :3p e S such that \\x - p||oo < e} 
i.e. the tube of radius e of S in the || • || norm. Then 

Pe{S) C T^{S) 

Proof. Suppose x G PeiS). Then 3p e S such that x and p he within the same e pixel. But by the 
definition of pixels, this implies that \\x — p\\oo < £■ Therefore x G T^{S). □ 



10 



BRANDON ROWEKAMP 



2. PiXELATIONS OF FUNCTIONS AND FIRST-ORDER APPROXIMATION 

In order to reach our overall goal of recovering sufficiently nice sets from their pixelations, we 
must start by examining simple sets. The graphs of functions on compact domains provide a good 
starting point. We will find that in this case many analogues of basic theorems from calculus appear 
in a pixelated version, and these facts will allow us to approximate the graphs of these functions. 

Before we proceed in our investigation we need to introduce a basic vocabulary that will be used 
throughout the paper. 

Definition 2.1. Fix e > and a bounded set 5 C M^. 

(1) A point a G M will be called e-generic if x G R \ eZ. For such a point a we denote by Ie{a) 
the interval of the form [ne, (n + l)^], n & Z that contains a. 

(2) For an interval [a,b] c M we define the vertical strip 

S[a,fe] := [a,b] X M 
For every A; G Z we denote by i>e,k the vertical strip 

[ke, {k + l)e] x M = ^ke,{k+i)e] 
For any £-generic point a G M we denote by 2>£{a) the strip 

/e(a) X M = k := [x/k\. 

(3) A column of PeiS) is the intersection of Ps{S) with a vertical strip Se_fc, fc G Z. The con- 
nected components of a column are called stacks. 

(4) For every e-generic a G M, we define the column of a pixelation Pe{S) over a to be the set 

C,{S,a) := §e{a)nPeiS). 

In other words, (^^(S', a) is the union of the pixels in Pe{S) which intersect the vertical line 
{x = a}. When S is the graph of a function /, we will use the notation C£{f, a) to denote 
the column over a of the pixelation Pe{f)- 

(5) We define the top, bottom and respectively height of a column to be the quantities 

T,{S,a) := sup{y : (a,y) G Ce{S,a)}, 
B,{S,a) :=M{y:{a,y)eC£{S,a)}, 

and respectively 

h£{S,a) :=^(r,(5,a)-S,(5,a)). 

□ 

Since stacks are constrained to be within the same column, a stack can usually be thought of as 
two pixels together with all the pixels that fall vertically in between them. 

This language allows us to state an analogue of the Intermediate Value Theorem. 

Theorem 2.2 (Pixelated Intermediate Value Theorem). If f : [a,b] ^ M. is a continuous function, 
then for every x G [a, b] \ eL, the column Ce{f, x) consists of exactly one stack. 

Proof. We argue by contradiction. Suppose that the column C£{f,x) has at least two stacks. This 
means that within the column over x there are two stacks with a gap of empty pixels in between them. 
This implies that there exist an interval [ce, de], not contained in the range of /, and real numbers 
xi,X2 G Ie{x) such that 

f{x2) <cs <de < f{xi). 
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Since / is continuous, this contradicts the Intermediate Value Theorem. Therefore the column over x 
has only one stack. □ 

Note that the Pixelated Intermediate Value Theorem immediately implies that a cycle cannot appear 
in the pixelation of a continuous function. Therefore the pathological behavior witnessed in Appendix 
[a] will not occur in the pixelations of functions, which is our first hint that pixelations of 
functions are the nicest type of pixelation to work with. 

This analogue of the Intermediate Value Theorem tells us that we can think of the pixelation of 
a function that assign a stack of pixels to each value in the domain of the function. The stacks 
themselves are worth investigating. To generate a stack in the pixelation, the function must attain 
values near the top and bottom of the stack within that column. This means that the average change 
of a function over a column must be related to the height of the column. This relation will be shown 
in the following two theorems. The first tells us that the height of a stack is bounded by the derivative 
of the function. 

Proposition 2.3. Let / : [a, 6] — >• M Z^e a Lipschitz continuous function. We denote by ||/'||oo the best 
Lipschitz constant, i.e., 

x^y 1-^ y\ 

Then for any e > and any x G [a, 6] \ eZ we have 

he{f,x) < ||/'||oo+2. 



Figure 7 . If the line on the left were any shorter it would generate the pixelation 
on the left. Thus, the minimum range of a function in a column is two less than the 
height of the column. 
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Proof. We argue by contradiction. This means that there exist e > and x G [a, h] \ eL such that 

> ll/lloo + l. (2.1) 



Set n:= [a;/ej G Z so that ne < x < (n + \)e. The inequahty (2.1 1 is equivalent to the condition 

r,(/,x)-i?,(/,x) >e(||/'||oo + 2). 

This implies that 

3c, d G (ne, (n + l)e) such that |/(c) - > 1 1/| |ooe, 

because the pixelation can be one pixel higher or lower than the graph of /; see Figure [7] We have 
thus found c and d such that 

|/(c)-/(d)| ^ ||/||oo£ 

\c-d\ > e -11-^11- 

But this contradicts the definition of 1 1 /' 1 1 oo ■ □ 

Remark 2.4. (a) Observe that if x G [a, h\ n eL then 

/i,(/,x) < 2(||/'|U + 1). 

(b) The above theorem can be used to bound the height of a stack by the derivative of / over that 
stack by considering the restriction of / to that column. □ 

The next theorem suggests that the height of a stack is close, in a crude way, to the absolute value 
of the derivative in the column for sufficiently small e. The intuition is that if / is C^, then any 
sufficiently small section of its graph will be very close to the graph of a line, and the height of a 
stack in the pixelation of the graph of a line must be very close to the slope of that line. 

Proposition 2.5. Let f he a function on [a, 6]. Then, for any e > and any xq G [a, 6] \ eL we 

have 

l/MI </i.(/,xo)| + ||/"||ooe + 2. 

Proof. Recall that /^(xo) indicates the interval of the form [ne, (n + n G Z, that contains xq. 
Let L be the linearization of / at xq, i.e., the function 

L:M^R, L(x) := (x-Xo)/'(xo)+/(xo). 

From Taylor's Theorem we deduce that 

llfll 1 
|i(x)-/(xo)| < ^^(x-xo)2 < -lirilooe' <e, Vx G /,(xo) n [a, 6]. 

This implies that 

r(L,xo) <r(/,xo) + ^|iri|ooe' + e, i?(L,xo) >S(/,xo)-^||/"||ooe'-e. 

Hence, 

/i,(L,xo) </i(/,xo) + ||/"||ooe + 2. 
On the other hand, L changes by /'(xo)e over the width of a pixel, which implies 

|/'(xo)| < he{L,Xo) 

□ 
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We would now like to approximate / through its pixelations. An intuitive way to do this is by 
creating a PL function whose graph lies within Pe{f), since this approximation would generate a 
similar pixelation. However, in practice this method requires making many choices, and it can be 
cumbersome to check that the graph of the function lies within the pixelation of /. A better approach 
is to form a PL function by connecting points within Ps{f) by straight line segments. In order to 
preserve higher order properties of the function, we must be careful when choosing these points. 

For example, suppose that the vertices of the PL approximation are centers of the top pixels in 
each column. Then the slopes of each line making up the PL function would have integer slopes, and 
could never approximate the derivative of / arbitrarily well. We could increase the accuracy of the 
approximation by connecting the maxima of every other column, or every third column. However, 
doing this decreases the speed in which the approximation converges point wise to /. This is the 
central tension involved in approximating a function from its pixelation. 

Before we further describe how to approximate a function from its pixelation, we show how to at 
least find the average change of the function, i.e., the slope of the secant line connecting the first and 
last points of /. 

Theorem 2.6. Let f be a function on the interval [a, 6]. Fix e > 0. Let c, d G M such that 

ia,c)GPe{f) and {b^d) € P,{f). 



Then 
and 



ge{x) := (x - a)^ — - + c. 



|/(x) - geix)\ < 7(||/'||oo + 2)e + 3||/'||oo(6 - a) + ||/"||oo(6 - «)' (2.2) 
If'ix) - g'Ax)\ = If'ix) -p^\<ib- a)\\f"\U + ^^"{1"°°!^^' , (2.3) 



for all X G [a, b]. 



Proof. Note that it suffices to prove the above inequalities in the special case when x G [a, b] \ eL. 
For ease of notation let 

M:=||/'|U, iV:=|iri|oo, ^ 

b — a 

and note that is the slope of the secant line to / over the interval [a, 6]. Define 

d — c 
b - a 

and note that this is the slope of g^- Note that c and d differ from /(a) and f{b) by at most the height 
of a stack. 

Lemma 2.7. With m f and irig defined as above, 

K-™,|<W1U±2)£^ (2.4, 



Proof of Lemma 2.7 Theorem |2 . 3 1 implies that a height of a stack in Pe{f) is at most M + 2 so 

\f(a)-c\,\f{b)-d\ < (M + 2)e. 

This immediately implies 



□ 
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The Mean Value Theorem implies that 3^ E [a, b] such that f'{C) = "^J- Let h{x) be the lineariza- 
tion of / at ^,1-6-, 

h{x) := (x-e)m + /(e). 

Thus we have 

\9e{x) - h{x)\ = \{x - a)mg + c-{x- - /(^)| 

< |(:c-a)mg-(x-0"i/| + |d-/(OI (2.5) 
First we will bound the first term on the right hand side of ( |2.5| ). 

|(x — a)mg — (x — ■^)?Tij| = |(x — a)mg — (x — £,)^g + {x — ~ ~ 

< |mg|(|x - a\ + \x - ^1) + \mg - m/|(|x - ^|) 
To proceed further we need to use the following facts. 

• |x — a|, |x — ^1 < |6 — a|. 

I I ^ 2(M+2)e 

• \mg-mf \ < \_J . 

• Since |mj| < M, we have 

mg < M + ^- 

b — a 

Therefore 

,/ N / , X.,. 2(M + 2)e, „ ,2(M + 2)e 
|(x - a)mg - (x - < 2(6 - a)(M + ^- ^) + (6 - a)^- ^ 

Qi (X 

= 2M(6- a) + 6(M + 2)e. 



Now we bound the other term from (2.5 1, |c — /(^)| = |c — h{(^)\. We have 

|c - /(Ol < |c - /(a)| + |/(a) - hia)\ + |/i(a) - 

Note that \f{a) — c\ < {M + 2)s, since stacks are at most M + 2 pixels tall. Taylor's Theorem implies 
that 

|/(„)-M.)l<^^. 

Finally \h{(,) — h{a)\ < M{b — a), since h has slope at most M, and |^ — a| < (6 — a). We conclude 
that 

1/(0 - c| < (M + 2)e + + M(6 - a). 

Combining all of these bounds in ( |2.5[ ) we deduce 

|5.(x) - h{x)\ < 3M(6 - a) + 7(M + 2)e + 

Finally, 

\f{x)-h{x)\ < \ , 

so that 

|/(x) - g{x)\ < 7{M + 2)e + 3M(6 - a) + iV(6 - af. 
which is the bound ( |2.2[ ). 

The bound ( |2.3| ) is obtained as follows. 

- Ssla;)! = |/'(x) - mg\ < \mg - m/| + \mf - f'{x)\ 



'i0 2iM^ 2(M^^^^ 

b — a b — a 
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□ 

This Theorem tells us that the slope of a secant line to / can be accurately approximated up to its 
derivative (by restricting / to various intervals we can approximate any secant line.) We would like 
approximate / through the use of multiple secant lines. However to accurately approximate a small 
secant line, we will need e to be very small. On the other hand, to accurately approximate / we will 
need to use many (and therefore small) secant lines. Thus we will need to know the size of a secant 
line that can be accurately approximated for a given e. Before we can describe how to go about this 
process, however, we will need additional terms. 

Definition 2.8. Fix e > and a bounded set S. 

(1) An e-profile of is a set of points in the plane with the following properties. 

(a) Each point in fig is the center of an e-pixel that intersects S. 

(b) Every column of Pe{S) contains precisely one point of n^. 

(2) A top/bottom e-profile is a profile consisting of the centers of the highest/lowest pixels in each 
column of Pe{S). 

(3) An £-sample of 6" is a subset of an e-profile. 

□ 

Definition 2.9. Suppose pi, . . . ,pjv is a finite sequence of points in M^. (The points need not be 
distinct). We denote by 

{pi,P2, ■ ■ ■,Pn) 

the PL curve defined as the union of the straight line segments \pi,p2], ■ ■ ■ , \pn-i,Pn]- D 

Observe that each e-profile Ilg of a set is equipped with a linear order ^. More precisely if pi,P2 = 
{x2, 2/2) are points in n^, then 

Pi ■< P2 <^ x{pi) < x{p2), 

where a; : — > M denotes the projection (x, y) i-)- x. In particular, this shows that any e-sample of 
S carries a natural total order. 

If E is an e-sample of S, then the PL-interpolation determined by sample S is the continuous, 
piecewise linear function L = L^. obtained as follows. 

• Arrange the points in E in increasing order, with respect to the above total order, 

^ = {^0 ^ 6 ^ 6 ^ • • • ^ Cn}, n + 1 = 

• The graph of is the PL-curve (^o, ^1, • • • , Cn)- 

In apphcations, the sample sets S will be chosen to satisfy certain regularity. 

Definition 2.10. (1) If cr is a positive integer and is an e-profile, then an e-sample with spread 
(7 is a subset 

S = {^0 ^ • • • ^ ^n} C U,{S) 

such that the following hold. 

• The points .^0 and ^„ are the left and rightmost points in the profile. (That is for each 
p e n^, x(Co) < xip) < x{Cn)-) 

• For any p G lie, there exists ^ G S such that \x{p) — x{$,)\ < ea. 

-a < -\x{^k) -x{Ck-i)\ <(^, VA; = l,...,n. 
z e 
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(2) A spread is an increasing function e i— )• a{£) from the positive real numbers to the positive 
integers with the following properties: 

(a) lim£\^oo-(e) = oo 

(b) lime\^oeo-(e) = 

□ 



The next theorem tells how to approximate / using profiles and spreads while maintaining first 
order properties. 

Theorem 2.11. Let f : [a,b] -^M.be a C^-function and a{e) be a spread function. Fix e > and let 
H be an e-sample of the graph of f with spread cr(e), 

S = {Co -< 6 -<•••-< Cn}, Ck = {xk,yk), k = 0,l,...,n. 
Denote by /= : [xq , x„] — )• M the PL-interpolation determined by S. Then 

I f{x) - fs{x) I < 7( ll/lloo + 2)e + 3\\f'\\^ea{e) + \\f"\\oo{ea{e)f , Vx G [a, b] n [xo,x„], 



and 



\f'{x) - /4(x)| < + + \\f"\Ua{e), 



for all X G [a, b] n [xq, Xn] \ {xq, • • • , Xn}- In particular, if 

limefjfe) = and limufe) = oo, 

then e \j 0, converges to f in the Sobolev norm W^'^ for all p G [1, oo). 

Proof. Let x G [xq, x„] and A; G {0, 1, . . . , (n — 1)} such that x G [xfc, x^+i]. On this subinterval, 
the function /= is defined by 

/h(x) := (x-Xfe) Wk- 

Xk^\ - Xk 



Thus /= is a function of the type described in Theorem 2.6 where 

a = Xk, b = Xk+i, c = yk, d = yk+i- 

Noting that 

^e(j(e) <b — a < e(T(e), 

we obtain the desired error bounds. 

Now suppose that linigx^o = and lim^x^o f C^) = c»- Then note that every term in both 

the error bounds of | / — /s| and | / — /^l is a constant multiplied by e, or ea{e) all of which 
converge to as e \ 0. This implies that 

rb 

lim / (|/(x)-/=(x)|P + |/'(x)-/4(x)|P)dx = 0, Vl<p<cx). 
Therefore /= converges to / in the Sobolev W^'^ norm. 

□ 



This theorem allows us to approximate functions using only their pixelations. Note that any 
profile can be chosen for the approximation. We can see this in the pixelation by noting that the top 
profile and bottom profile both converge to the graph of the function, so any profile chosen in between 
will be squeezed onto the graph. 
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To ensure the correct convergence on the derivative, we only need the two basic features of a 
spread: that as e \ 0, a{e) — > oo and ea{e) \ 0. In general we would like the spread to increase 

quickly, so as to get an accurate approximation on the derivative. The two conditions on the spread 

-1 

mean that the quickest increasing spread will look like [e 1+^^ ] where r > 0. 
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3. Approximations of Elementary Regions and Curvature Approximation 

In this section we use the properties of pixelations of functions to study the pixelations of simple 
two dimensional sets. We deal with the simplest case: that of a region between the graphs of two 
functions. 

Definition 3.1. A subset 5" C is said to be elementary (with respect to the x-axis) if its can be 
defined as 

S = 5(/3,r) := :xG [a,6],/3(x) <x<r(x)}, 

where /3, r : [a, 5] — > M are functions such that /3(x) < t{x), Vx € [a, h]. The function /3 is called 
the bottom of S while r is called the top of 5. Note that this includes the situation where /3 = r, in 
which case Pe{S) = Peif3) = ^^(r). □ 

In the remainder of this section S will indicate an elementary set. We first note that like the 
pixelation of a function, each column of Pe{S) contains only one stack: 

Proposition 3.2. If S = S{(3,t) is an elementary set, then for every x G [a, 6] \ eZ, the column 
C£{S, x) consists of exactly one stack. 



Proof. Fix an e-generic x G [a-,h]. By Theorem 2.2 the columns Ce{f3,x) and Cs{t,x) consist of 



single stacks. If these two columns intersect, then the conclusion is obvious. If they do not intersect, 
then any pixel in the strip Ss{x) situated below the stack Cs{t, x) and above the stack Ce(/3, x) is a 
pixel of Pe(S'). This again proves that the column Ce(S', x) consists of a single stack. □ 

Definition 3.3. Fix e > and an elementary set S = S{(3, r). 

(1) An s-upper/lower profile of 5 is a profile of the e-pixelation of the top/bottom function. An 
e-upper/lower sample is a sample of an upper/lower profile. 

(2) An e-upper profile 11+ is said to be compatible with an e-lower profile 11^ if p+ G 11+ and 
p~ G n~ lie in the same column of Pe{S) we have p+ G 11+ and p" G 11^ lie in the same 
column of Pi; (S) we have 

(3) An e-upper sample H+ is said to be compatible with an e-lower sample E~, if the following 
hold. 

• There exist compatible e-upper/lower profiles 11+ such that H+ C 11+ . 

• A strip §£ fc contains a point in H+ if and only if it also contains a point in E~ . 

(4) Suppose that H+ are compatible upper/lower samples of S 

The PL-approximation of S determined by these two samples is the PL-set bounded by the 
simple closed PL-curve 

We denote this approximation by | S{Ej , H+ ) | . □ 

These definitions are extensions of the ideas of profiles, samples and PL approximations found 
in section 2. Note that the upper and lower profiles are defined from the pixelations of the top and 
bottom function, which might not be known. However, from the definition of an elementary set, every 
pixel on the top of a column must intersect the pixelation of the upper function and similarly every 
pixel at the bottom of a column must intersect the lower function. Therefore, an upper profile can be 
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Figure 8. An example of PL approximation. The left figure shows a pixelation of 
an elementary subset, together with a compatible upper and lower profile. The upper 
profile is indicated in blue while the lower profile is in red. The right figure shows 
the PL approximation of compatible samples of these approximations with spread 
2. The PL approximation for the combined upper and lower samples would be the 
region in between the upper and lower approximations. 

thought of as a profile chosen from the top pixels of Ps{S) and a lower profile can be thought of as a 
profile chosen from the lower pixels of Pe(<S'). 

In section 2 we noted that the PL approximation preserves first order properties for a good choice 
of spread. For elementary sets, we would like to prove something stronger. Namely, we would like to 
show that the total curvature of the set is preserved for good choice of spread. 

The total curvature of a immersion 7 : [a, b] — )• M? is defined as follows. Define a C^-function 
■.[a,b] such that 

^7(i) = (cos(0(i),sin(0(i))), 
where a dot ' denotes the t-derivative. We set 

<t) ■■= m\. 

The scalar K{t) is called the curvature of 7 at the point ^{t). We define the total curvature of 7 to be 

K{j) := [ K{t) dt. 



Suppose now that 7 : [a, b] — )■ continuous and piecewise C^-immersion, i.e., there exist a finite 
subset {to, . . . , tiy} C [a, b], ij^i < U, for 

a = to < *i < • • • < tu-i <ty = b 

such that the restriction 7, := 7|[[j._j^f.] is a -immersion for any i = 1, . . . ,v. The curvature of 7 
at a jump point 7(ti) is the quantity 

^^{U) = \e{tt) - e{tT)\ := I lim e{t) - llrn e{t)\. 
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We define the total curvature of 7 to be 



i=l i=l 

fo, 7(&)/7(a) 
\\ei4)-9it~)\, 7W=7(a)- 

For more details we refer to | [T8l and ll20l §2.2]. 

Proposition 3.4. Let h : [a,b] ^ be a function. Fix a spread a satisfying the properties 

lim e(j(e)^ = 00, (3.1a) 

lim ecT(e) = 0. (3.1b) 
For every e > let be a e-sample spread a. We denote by be the PL interpolation of S^. Then 

\unK{Th,) = K{Th), 

where for any function f we denote byTf its graph. 

Proof Let be the ordered set of x-values of the sample E^, 

= {xo < xi < ■ ■ ■ < Xn} C [a, b], 

and set 

Pi := {xi,h{xi)) e Th. 
Define to be the PL function with graph {pQ,pi, . . . ,Pn)- We will prove two things. 

lim{K{he)-K{se)) = 0, (3.2a) 

e-s>0 

lim K(se) = K(h). (3.2b) 

£-5>0 



Proof of (3.2a I. For 1 < z < n let indicate the slope of the i-th line segment of he and mf 
indicate the slope of the i-th line segment of Sg, i.e., the slope over the interval Xj]. 

Since has spread a, the width of the interval Xj] is at most ea. Then since is the slope 
of a line connecting two points in the sample H^, and m| is the slope of the secant line over the same 
interval, Lemma [277] implies that 

\ h Sl^r,ll^ll°°~'~^ /TON 

|m- - I < 2 — (3.3) 

In particular, \m^ — mf\ — as e — )■ since cr — t- 00. 

Define 9^ G (— |, |) as the angle that the {i + l)-st line component of makes with the x-axis 
and likewise 9f G (— f ) f ) as the angle that the {i + l)-st line component of s makes with the x-axis. 
More precisely, 

tan(6'f ) = rui, tan(6'|) = m-. 
Therefore, the difference formula for tangent implies: 

h s 



1 + m'^mf 

so that 



tan(6l- 



— m; 



arctan ( — ^ ) . (3.4) 
1 + niimf 
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Now, the estimate (3.3 1 implies that for small e, we have 



my 



1 + m'^mf 



< 4- 



\h'\\oo + 2) 
a{e) ■ 



Indeed, the denominator is at least ^ because, for e sufficiently small, we have 

1 + = 1 + ( m,'^ + O {a{e)-^) ) 

= 1 + {m^f + m^O {a{e)-') = 1 + {m'^f + O {\\h'\\^aie)-'] 
Then, using the Taylor expansion of arctan around 0, we have 



Of = o 



1 



(Tie 



(3.5) 



Recall that 



and we deduce that 



i=l 



i=l 



\K{Se)-K{h,)\<Y,m + \5i- 



Note that since has spread a, 
and (|3.5|) implies 



n = 



1=1 
1 

ecr 



tm + \s,-A)-o[^). 



Since we have assumed that linigx^ ea{e)'^ = oo, the last estimate implies (3.2a I. 

/ the equation: 

[l,h\x)) = (cos(0(x)),sin(e(x))) 



Proof of (3.2b I. Define the function 6 by the equation: 

1 



\l + {h'{x)f 

The function 6{x) gives the angle between the x-axis and the tangent line of h at the point x. As 
before let m| indicate the slope of the i-th line segment of and Of be the angle that the f-th line 
segment of makes with the x-axis. The Mean Value Theorem implies that for each i there exists a 
on Xj] such that 

ml = h\i.i) 

It then follows that 

et = o{i^) 

since tan(m|) = Of and tan(/i'(^j)) = 0{^i). Furthermore, for each i > 1, the Mean Value Theorem 
guarantees the existence of a Ci on ^j] such that 

Recall that the total curvature of Se is calculated by the summation: 



i=2 
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Using the definitions of and Q we find tliat 



K{ss) = j2\m)-o{^i-i)\ = Y. 



1=2 



i=2 



iii-ii-l) = Y.\^'{C^m^-^^-l) (3.6) 



i=2 



Note tliat last sum in (3.6 1 is almost a Riemann sum for \6'\. It is not a Riemann sum because the 
values {■^1,^2, • • • , in] do not form a partition of [a, b] (since > a and ^„ < h.) However, it is true 
that: 



\\m[\d'{a)m-a) 

e\0 \ 



+Y.\^'{Ci)m-ir-i) + \o'm{h-in)] = / \e'\ 

i=2 J 



x)\dx 



(3.7) 



Both (6 — in) and (^i — a) are intervals of size at most ea{e) (since our sample was chosen with 
spread a.) Therefore both of these terms vanish as e — )• (by our constraint on a.) Therefore equation 
(|3.7|) simplifies to 



[x)\dx 



(3.8) 



Combining (3.8 1 with (3.6 1 shows that 



lim Kisr 

e\0 ^ 



r{x)\dx = K{h) 



which proves (3.2b i. 

Given both (3.2a) and (3.2b I we simply note that as e — )• approaches K{s£) which in turn 

approaches K{h). □ 

In applications, it will be necessary to deal with sets which resemble elementary sets, but whose 
boundaries are possibly only piecewise C^. Let us define precisely the notion of piecewise C"^ func- 
tion. 

Definition 3.5. A function / : [a, 6] — )■ M is said to be piecewise if there exists a finite set 

S = {a = so<si<---<s^ = 6}, 

such that the following hold. 

(1) The function / is continuous. 

(2) For any j = 1, . . . ,i, and any k = 1,2 the restriction of / to the open interval (sj-i, Sj) is a 



C"^ function and the limits 



lim f^^\x), lim 



exist and are finite. 



We say that S is the singular set of / and that the points Cj = {sj, f{sj)), j = 0, ... ,^ are the 
corners of the graph of /. The integer £ is called the length of the singularity set. □ 

Proposition 3.6. Suppose that f : [a, 6] — ?■ M is a piecewise function with singular set 

S := {a = So < si < ■ ■ ■ < S£ = b Y 

Fix a spread satisfying the properties l \3.1a I and (S.lbl For every e > let be an e-sample with 
spread a satisfying {3.1a) and {3.1b). For any e > we denote by fs be the PL interpolation ofE^. 
Then 

\m.K{Tf^) = K{Tf). 

£\0 
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Proof. We will follow the same strategy as in the proof of Proposition |3 .41 with some expected mod- 
ifications due to the presence of singularities. For the clarity of exposition we will assume that the 
length £ = 2 so that 

S = { a,b,si; si £ (a, 6) }. 

Set 

/ := f\[a,si], '■= f\[si,b]- 

Let 

= {{xj,yj); j = 0,...,n = n{e)}. 
We denote by the ordered set of x-values of the sample H^, 

Xs = {xQ < xi < ■ ■ ■ < Xn, n = n{e)} C [a, b], 

and set 

Pi := {xi,f{xi)) £ Th- 
Define g,, to be the PL function with graph {pQ,pi, . . . ,pn). We will prove two things. 



lhn(K(f,)-K{g,)) = 0, 
e-!>0 



\imK{g,)=K{f). 



(3.9a) 
(3.9b) 



Proof of (|3.9a|l. For 1 < i < n let mf denote the slope of the i-th line segment of and denote 
the slope of the i-th line segment of (7e, i.e., the slopes of the segments over the interval Xj]. 
Since has spread a, the width of the interval is at most ea. 

Since the function / is Lipschitz continuous we deduce from Theorem 2.3 that 

\f{xi)-yj\ < (||/'||oo + 2)e. 

This implies that 

I f 01 ^ n 1 1^ 1 ^ 

m- — < 2 

In particular, \m{ — mf\ — )• as e — )• since cr(e) oo as e 0. We can now conclude as in the 
proof of ( |3.2a i. 

Proof of (3.9b I. For every e > there exists a unique is G {!,..., n{e)} such that 

Xi^-l < Sl< Xi^. 

We set = /(si), p* := (si, y*), and we denote by hg the PL-function with graph 

{PO, ■ ■ ■ ,Pi,-l,P*,Pi,, ■ ■ -Pnie))- 

We set 

hg ■= ^e|[a,si]; := /lelfs^^^fo], 

We denote hy 9^ G [0, vr) the angle between the two nontrivial line segments of the graph of h that 
have as common vertex. Then 

i^(r,j = i^(r,-) + i^(r,+) + e,. 



From Proposition 3.4 we deduce that 



iimK(r,±) = i^(r.±), 
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while 9e converges as e — )■ to the angle between the left and right tangents to the graph of / at , 
so that 



iimi^(r,j = i^(r^). 



We thus have to prove that 



iini(/^(r,j-E:(r,j) =0. O.io) 

To analyze the difference K{Tg^) — KiTh^) we distinguish two cases. 

A. si = Xi^. In this case = he and thus K^Tg^) — KiTh^) = 0. 

B. si < Xi, i = is- We introduce the following notation from which we suppress the e-dependence. 

• a_ G [0, tt) denotes the angle between the vectors pj„2Pi-i Pi-ip\- 

• a+ G [0, tt) denotes the angle between the vectors pi-ip\ and piPi^l. 

• /3_ G [0, vr) denotes the angle between the vectors pi-2Pi-i and pi-ipl. 

• /3+ G [0, vr) denotes the angle between the vectors pTp^ and pjpj+i. 

• ^£ G [0, vr) denotes the angle between the vectors Pi-ipl and p^:p\. 

• denotes slope of the line pi^ip^. 

• m+ denotes the slope of the line p^pi. 

• = denotes the slope of the line pi^ipi. 

Then 

K{Tg^) - K{Th^) = (a_ + a+) - + /3+ + 6^), 



arctanm^_j_]^ — arctanm^ 



a± = |arctanmf_j_]^ — arctanm*| , /3± 
Now let us observe that 

mf^i =mt + o(l), 

where, following Landau's convention, we denote by o(l) a quantity that goes to zero as e — )• 0. Thus 

/3± = o(l), 

We deduce that 

K{TgJ - K{ThJ = a+ + a.-ee + o(l) 

= |arctanm~ — arctanm*| + [arctanm^ — arctanm,,,! — + o(l). 

We now remark thafl 

larctanm" — arctanm^,! + larctanm^ — arctanm*! = (3.12) 



The equality (3.12 1 is a classical fact of Euclidean geometry. More precisely, in the triangle Pi-ip*Pi, 
the sum of the interior angles at the vertices and pj+i is equal to the sum of exterior angle at p^; 
see Figure |9] Using ( |3.12 ) in ( |3.11| ) we obtain (3.10 1. This completes the proof of (3.9b I, and thus 



the proof of the Proposition |X6] □ 

This will allow us to handle the case of piecewise linear sets. Therefore the following corollary 
will be useful. 



^Compare with Lemma fTSi Lemma LI]. A word of warning: while the main conclusion of that Lemma is true (total 
curvature never decreases upon adjoining a point to to a PL-curve, the claim that it does not change for planar curves is 
not true. 
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Pi-i 




x+y+z=7r 

9e+Z=7t 



Figure 9. A classical fact of Euclidean geometry 

Corollary 3.7. Let S{I3, r) be a set which can be expressed as 

S = {{x,y) : X ^ [a, 6], /3(x) <y < t{x) ] 
where /3, r : [a, 6] — )■ M are piecewise functions such that 

f3{x) < t{x)^ Vx G (a, 6). 



Fix a spread a{e) satisfying the conditions {3.1a) and {3.1b). For each e we choose compatible e- 
upperAower profiles "Ef of S with common spread (j{e). We denote by the PL approximation of 
S defined by these samples. Then 

\\uiK{dSe) = K{dS). 

e\0 



Proof. Let 

The compatibility condition implies that 



^ }) ^"k — i^^^Vk^ 



^k =H=:xk, < y^, V/c = 0, 1, . . . , n. 

Let be the PL function whose graph is the bottom part of the boundary of Ss, and to be the PL 
function whose graph is the top part of the boundary of S^. Let (e) indicate the slope of the i-th 
line segment of the graph of Pe and similarly let mj{e) indicate the slope of the i-th line segment of 
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the graph of r^. We have 

K{dSe) = Ivr — arctan(m^(e) 



+ I arctan(mf (e) ) — arctan(mf_-|^(e) )| 

i=2 

+ I arctan(m(^(e)) — 7r| + |7r — arctan(my)| 

n 

+ I arctan(mj'(e) ) — arctan(m[_]^(e) )| 

1=2 

+ I arctan(m^(e) ) — 7r| 

which can be rewritten as 

K{dSe) = Ivr — arctan(m^(e) )| + | arctan(m^(e) ) — 7r| 

+ Ivr — arctan(m[(e) )| + | arctan(m^(e) ) — vr| (3.13) 

+ K{Pe) + K{Te) 

Proposition |3.6| implies 

\\mK{Pe) = K{P) and lim i^(re) = K{t). (3.14) 

e\0 e\0 

Now note that each line segment is defined by connecting two points in the pixelation of /3 or r over 
an interval of at most ecj(e). Since as e — 0, ecr(e) — 0, Theorem |2. 6| implies that 



lim m^(e) = /3'(a), lim m^(e) = 
limm[(e) = /(a), lim?Ti^(e) = 

£\0 e\0 



(3.15) 



Combining ( [37T3] ), pi?] ), ( [STTS] ) we find that 

lim K{dSe) = |vr — arctan(/3'(a))| + | arctan(/3'(6)) — 7r| 

e\0 

+ |vr — arctan(r'(a))| + | arctan(r'(6)) — vr| (3.16) 
+ K(/3) + K(r). 

Note that |vr — arctan(/3'(a))| is the value of the angle between the vertical line x = a and the tangent 
line of j3 at a. Similarly each other difference on the right hand side of ( 3. 16 1 corresponds to an angle 
at one of the corners of dS. Therefore the right hand side of the (3. 16 1 is equal to the K{dS), so the 
corollary holds. □ 
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4. Approximations of PL Sets and Morse Theory for Pixel ations 

In previous sections we have dealt only with the simple regions and we investigated mainly geo- 
metric properties of these regions and their pixelation. In this section we turn our attention to the 
relationship between the topology of a PL-set and those of its pixelations. 

Surprisingly, this is a nontrivial matter. In general, the homotopy type of a region may not be 
not preserved when taking the pixelation. Worse, the homotopy type may never be recovered in any 
pixelation, even for small resolutions s. This type of bad behavior can happen even for a simple 
PL case. Consider the set S composed of the rays starting from the origin and proceeding in the 



positive direction with slopes ^ and ^ (see Figure 10 1. A careful examination of Pe{S) reveals the 



existence of cycle. Worse yet, the pixelations for smaller values of e are simply contractions of the 
larger pixelations. This means that Pe{S) contains a cycle for all small e. For a taste of how much 
worse can this get we refer to Appendix [A} 









- i 







Figure 10. The pixelation of the lines with slope ^ and j. After this point the 
pixelations of the upper and lower lines diverge permanently and no more holes are 
formed. 

First, some good news. The next two results imply that if two compact planar sets are disjoint, then 
their pixelations will also be disjoint for sufficiently fine resolutions. 

Theorem 4.1. Let S cM? be a compact set and let x G such that x ^ S. Set 

m := inf d{x, s), 

where d{», •) is the Euclidean metric on R^, so that m is the Hausdorff distance from x to S. Then 

Ve < f , X ^ Pe{S). 

Proof. For x to lie in Pe{S), both x and a point of S must lie in the same pixel. However the furthest 
two points can be apart in a square of side length e is (if they lie on opposite corners of the 
square.) But since e < ^, ey/2 < ^^^^ < m, so it is impossible for a point of S to lie in the same 
pixel as x. Therefore x ^ Ps{S). □ 

Corollary 4.2. Let S cM? be a compact set with finitely many connected components. Then for suf- 
ficiently small £, the number of connected components ofPs{S) agrees with the number of connected 
components of S. 
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Proof. First note that Pe{S) cannot have more connected components than S, since Pe{S) contains 
only pixels that intersect S. 

Since S is compact with finitely many components, there is a minimum positive distance 5 between 
components. For each component K of S define 

Tk ■■= {x G : inf d{x, k) < -} 
fcg/< 2 

and note that its topological boundary is 

dtopTK = {x G : inf d{x, k) = -} 
kaK I 

Then each point in dTx is | away from K in the Hausdorff metric. Furthermore, since 5 is the 
minimum distance between components of S, each point in BTk is at least | away from any connected 
component. Therefore, for e < g we have for each K 

dTK n Pe{S) = 

This implies that for sufficiently small e each connected component of Pe{S) intersects at most one 
connected component of S. However, since 

S C Pe{S) 

it is also true that each connected component of S is contained within one connected component of 
Pe{S). Therefore for sufficiently small e, Pe{S) and S have the same number of connected compo- 
nents. □ 



The corollary guarantees that the zeroth Betti number of a set coincides with those of its sufficiently 



fine pixelations. Theorem 4.1 also suggests that, for small e, the only way that the topological type 
of Pe{S) can disagree with that of S is if Pe{S) has additional cycles (since every "real" cycle will 
show up in the pixelation). We will refer to these cycles as holes. More formally the holes are cycles 



in Pe{S) that are not contained in the image of the inclusion induced morphismlations. Theorem 4.1 
also suggests that, for small e, the only way that the topological type of Pe{S) can disagree with that 
of S is if Pe{S) has additional cycles (since every "real" cycle will show up in the pixelation). We 
will refer to these cycles as holes. More formally the holes are cycles in PeiS) that are not contained 
in the image of the inclusion induced morphism 

Hi{S,Z)^Hi{Pe{S),Z). 

Thus recovery of S from Pe{S) will depend on distinguishing the cycles of PeiS) that correspond to 
real cycles from 5 from those that are merely artifacts of the pixelation. To discard these holes, we 
adopt a strategy inspired from Morse theory. 

Definition 4.3. Let S* C be a compact set and e > 

(1) For every e-generic x we set 

ne{x) = ns^s := # of connected components of Ce{S, x). 

(When the set S is understood from context we use the simpler notation instead of £.) 
We will refer to ne(x) as the stack counter function of S. 

(2) We define 

n(x) = ns{x) := # of components of { y G R; (x, y) G 5 }. 
We will refer to ^^(x) the component counter of S. 
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(3) A jumping point of ns is a real number xq such that 

n(xo) 7^ n{xQ ) := lim n(x) or n(x) / "-(a^o^) := lim n{x). 

We denote by the set of jumping points of ns- We will refer to as the. jumping set of 5. 

(4) A jumping point of is a real number xq G such that 

ns(xo) / n(x^). 

We denote by ds,e the set of jumping points of ns^e- We will refer to it as the e-jumping set 
ofS.' 

□ 

The function tells us how many stacks are in a column. A cycle has two "walls" and a "gap." 
That is to say, if there is a hole in ^£(5"), there will be a set of columns which have more components 
than their neighbors, since for a hole to close up stacks must overlap. Thus jumps of are a first 
indicator of the presence of cycles in Pe(<S'). To decide whether they are holes, as opposed to cycles 
coming from S we will rely on our next key technical result which is a substantial refinement of 
Theorem 14. II 



Theorem 4.4 (Separation Theorem). Let f, (7 : [a, 6] — )■ M be two Lipschitz continuous functions such 
that f{x) < g{x), Vx G [a, h]. Fix e > and suppose that 3xo S [a, 6] such that 

g{xo) - f{xo) < g{x) - f{x), Vx G [a, b], (4.1) 

Then for any e > such that 

( 3 + min(||/'||oo, ll^lloo) )e < gixo) - /(xq) (4.2) 

and any e-generic x € [a, b] \ eTL the column Ce{G, x) has two components. In other words, if 

mill ( g{x) — f{x 

x&[a,b] 

then for any e-generic x G [a, b] we have 



mill (^(x) -/(x)) > (3 + min(||/'||oo, ||fi''l|oo))e, 



nG,e(x) = nG(x). 



Proof Note that 

C,(G,x) = C,(/,x)uC,(5,x), 
and furthermore. Theorem |2.2| implies that each of these columns is connected. Therefore C^{G, x) 
will have two components only if C£{f,x) and C£{g,x) do not intersect. Since f < g and x is 
e-generic, this will occur when 

T,{f,x) <Be{g,x) 

or equivalently, 

B,{g,x)-Te{f,x)>e 

Now fix X G [a, b] and let i G Z such that ie < x < {i + l)e. Choose Xf,Xg G [ie, {i + such 
that 

/(x/) = max /(x), g(xg) = min o(x). 

Therefore we have 

B,{g,x)>g{xg)-e, Te{f,x)<f{xf) + e, 

so that 

Bs{g,x) -Teif,x) > gixg) - fixf) -2e. 
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We distinguish two cases. 

Case 1. Il^'lloo < ||/'||oo- We have 

Be{g, x) - Te{f, x) = g{xg) - g{xf) + g{xf) - f{xf) - 2e 

9{xg) - g{xf) + g{xo) - f{xo) - 2e 



> gixo) - f{xo) - e( 2 + H^'IU ) ^ e. 
Case 2. ||/'||oo < Ib'lloo- We have 

Beig, x) - T,{f, x) = g{xg) - f{xg) + f{xg) - f{xj) - 2e 

g{xQ) - f{xo) + f{xg) - f{xf) - 2e 



>ff(xo)-/(xo)-e(2 + ||/'|U) ^ e 



□ 



To proceed further we need to introduce some basic terminology. 

We define a convex polygon to be a compact set in that is the intersection of finitely many closed 
half-planes. Note that points, and straight line segments are examples of convex polygons. A PL set 
in M? is a finite union of convex polygons. Note that the topological boundary of a PL set is a finite 
union of straight line segments and points. A vertex of a PL set 5 is a point p on the topological 
boundary dtopS of S such that, for all r > sufficiently small, the intersection of dtopS with the 
closed ball of radius r and center p is not a straight line segment. We denote by Vs the set of vertices 
of a PL-set 5. 

Definition 4.5. (1) A convex polygon is a compact subset of M? which is the intersection of 
finitely many closed half-planes (note that line segments and points are examples of convex 
polygons). 

(2) A PL set (or piecewise linear set) in is a finite union of convex polygons. 

(3) A vertex of a PL set 5 is a point p on the topological boundary dtopS such that for all 
sufficiently small r > 0, the intersection of dtopS with the closed ball of radius r and center 
p is not a straight line segments. 

(4) For a PL set C M^, the set of vertices eVs is the collection of all vertices in S. 

(5) A PL subset in 5 C is called generic if for any two vertices pi,P2 G "^s 

x{pi) / X{p2). 

□ 



We will restrict our approximation technique to the PL case in order to simplify the conclusion of 
Theorem 4.4 Since our technique for determining holes will be motivated by applying Morse theory 
to projection onto the x-axis, we will also require the PL set to be generic (to avoid complications 
arising from clusters of critical points sharing the same critical value). Note that for a PL set S, the 
set of jumping points is contained in the set of x-coordinates of the vertices of S, that is 

as C {x{p); p^Vs). 

This inclusion could be strict. Take for example a V-shaped set, with the bottom vertex of the letter 
V situated at the origin. Then is not a jumping point of n. 
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Theorem 4.6. Let S be a generic PL set with jumping set ^g. Then there exist u = ^{S) G Z>o and 
^0 = £0(5') > 0, depending only on S, such that, ifO<£<eQ and x is an e -generic value such that 

dist(x, ds) > 

then ns^ei^) = i^s^x). 

Proof. Let xq < xi < • • • < be the jumping points of n = ns- We set 

Axi := Xi — Xi-i, Wi = 1, . . . , £, A := min Axj. 

i<i<e 

Note that n(x) is constant on each of the intervals (xi_i, Xj). For i = 1, . . . , £ we set 

Si := {{x,y) £ S; x G [xi_i, Xj] }. 
The set Si is a disjoint union of elementary sets (see Definition [3TT] for notations) 

S{l3i,j,Ti,j), j = 0,...,pi, 

"stacked one above the other", i.e., 

(3ifl{x) < Tifl{x) < A,l(x) < Ti,i(x) < • • • < I3i^p^{x) < Ti^p^{x), Vx G (Xi_l,Xi). (4.3) 



min nj. (4.4) 



From Proposition 3.2 we deduce that for any e-generic x G (xj_i, Xj) we have n(x) = pi. 

Each of the functions f3ij and Tj j is piecewise linear. Let rj j be the smallest widthj^ of a line 
segment of the graphs /3j j and Tij, and we set 

1 

Tj := - 

3 o<i<p, 

By definition, < ^Axj. We set 

5[xi_i + ri,Xi - ri] := { (x, y) G 5; x G [xi_i + ri,Xi - r^] }. 

The set ^[xi-i + rj_i , Xj — r^] is a collection of pi elementary sets stacked one above the other which 
have positive Hausdorff distance between them. Theorem 4. 1 implies that there exists 6i > such 
that, for e G (0, 5i) the pixelation ( ^[xj-i + r, Xj — r^] ) has exactly pi components. Therefore we 
have proven that 

35i > such that; Ve G (0, 6i) x G [xj_i + r^, Xj — rj] : n(x) = ne(x). (4.5) 

Set 

7/2 := Xi—i -\- Tij Zi Xi Ti. 

On the interval [xj_i, y^] each of the functions /3jj and Tij is linear and we denote by m~ (ftj ) and 
respectively m^(rj j) their slopes. For each x G (xi_i, Xj) and each j = 1, . . . ,pi we define the gaps 

7ij(x) = /3jj(x) - Tij_i(x), 7i.j := min{7ij(xj_i), 7j(yi)} = min 7ij(x), 

xelxi^i,yi] 

Fij := max{7ij(xi_i),7jj(yj)} = max 7i,j(x). 

xe[xi-i,yi] 

We plan to invoke the Separation Theorem |4.4| We want to prove that there exists Eq > and u > 
such that for e < we have 

min{7jj(x); Xj_i + ue < x < yi - ue} > (^3 + min( |m~(/3ij)|, |m~(rjj_i)| ) ^e. (4.6) 

Note that if 2z^e < Axj, then 

min{7jj(x); Xj_i + z^e < x < yj - i/e } = 7jj + z/|m"(/3ij) - m~(rjj_i)|e. 



The width of a line segment is the length of its projection on the x-axis. 
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e. 



(4.7) 



We can now rewrite (4^1 as 

> (^3 + min( |m"(/3ij)|, |m~(Tij_i)| )-z^|m~(/3ij) - m"(rij_i)| ^( 
To solve the last inequality we distinguish two cases. 

Case 1. = 0. In this case the slope of Pij must be different from the slope of Tij-i and we 
choose u = u^- large enough so that 

3 + min( |m"(^ij)|, |m"(Tjj_i)| )-i/|m"(/3jj) - m"(rjj_i)| < 0, 

e.g., 

3 + min( |m~(/3ij)|, |m~(Tij_i)| ) 



We then choose e~(i, j) small enough such that 2fj jeo < Axj, e.g., 

A 



+ 1. 



10m 



Case 2. 7^ ,,■ 7^ 0. In this case we choose eq = eo{i,j) small enough such that 



7iJ > (3 + min( |m {f3ij)\,\m (t;,j_i)| ) ^eq, 



e.g. 



^0 ihj) 



2(3 + min( \m |m(Tij_i)| ) 



Next we choose u = u- ■ such that 2veQ < Axj, e.g., 



A 



Finally, we define set 



£■ := min e v- = max v- ■. 

0<j<P^ ' 0<j<P^ ''-^ 



Theorem 4.4 implies 

3e~ > 0, z/r > such that Ve < e^, Vx G + i/^^e, Xi_i + rj] : n(x) = ne(x). (4.8) 
Arguing in a similar fashion we deduce 

3ef > 0, > such that Ve < ef, Mx G [xj — ri,Xi — f^^e] : n(x) = ne(x). (4.9) 
Now set 

Vi := maxfz^", u^), v := max fj, 

* * i<i<e 



Ei := min(e. ,ef,6i), Eq = min ej. 



Theorem 4.6 now follows from (4.5 1, (HTSh and (4.9 1 



□ 



This theorem tells us that jumps of occur within pixels from the jumps in n. A priori, it 
could be possible that, given a jumping point xq of n, there is no jump in within v{S) pixels of 

Xq. 



Theorem 4.7. Let S be a generic PL set, and Eq = £:o{S), v = 1^(5*) as in Theorem 4.6 Then, there 
exists El = Ei{S) such that if e < min(eo, ei) cind xq is a jumping point ofn = ns, then = ns^e 
has at least one jumping point in the interval [xq — ve, xq + ve]. 
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Proof. Since xq is a jumping point of n we have 

n(X([) / n(xo) or n(xo) / ri(xo ). 
We distinguish several cases cases. 

Case 1. '"-(xq ) > n(xo). Since 5 is compact, n(xo) > so that n(xQ ) > 2. For this to happen the 
vertical line oi x = xq must contain at least one vertex of S. Since S is generic, this line contains 
precisely one vertex of S, which we denote by po- 

We can find 5 > Q sufficiently small such that the interval [xq — 5, xq] will contain no new jumping 
points of 5. The set 

S[xo-5,xq] := G 5; x G [xq - 6,Xo] ] 

is disjoint union of simple types regions 

S{/3j,Tj), i = 0, . . . ,m = n(xQ ) - 1, 

"stacked one above the other", i.e., 

(3o{x) < To(x) < /3i(x) < ri(x) < • • • < /3m(x) < Tm{x), Vx G (xo - S, Xq), 

where l3j,Tj are piecewise linear functions. Since n(xo) < n{xQ) we deduce that there exists 
jo = I, . . . ,m such that 

^jo(^o) = Tjg-i{xo) and 7^ := /3j{xo) - Tj_i(xo) > 0, Vj / jo 

In particular, for any e > 0, the e-stack of S{PjQ,Tjg) over xq touches the stack of S'(/3j(,_i, Tj(,_i) 
over Xq. 

Now choose si sufficiently small so that for j / jo and e < ei, the e-stack of S{fij,Tj) over xq is 
disjoint form the e-stack of rj_i) over xq. Fixe < min(eo,ei). The above discussion shows 

that 

n(xo) = ne{xQ). 

Theorem 4.6 now implies that 

ne((xo - veY) = n(xo - ve) > n(xo). 
This proves that the interval [xq — re, xq] contains a jumping point of n^. 
Case 2. u^Xq) > n(xo). This situation can be reduced to the previous case via the reflection 

3 (x,y) ^ {-x,y) G 

Case 3. «.(xq ) < n(xo). The vertical line x = xq contains a unique vertex po of S. Moreover, this 
vertex has the property that there exists a tiny disk D centered at po such that the intersection of D 
with the open half -plane {x < xq} C is empty. In particular, this shows that po is an isolated 
point of the set 

Sx<xo = {i^,y) ^ S; x<xo}. 

If n(xQ ) = 0, the conclusion is obvious. We assume that n(xQ ) > 0. Choose 6 > such that the 
interval [xq — S,xo) contains no jumping point of S. Set 

R := S[x^_s^xo] \ {Po}- 

Then i? is a union of simple regions 

S{l3Tj), j = 0,1 . . . ,m = n(xQ ) - 1, 
where /3j and r, are piecewise linear functions such that 

l^oix) < ro(x) < /3i(x) < ri(x) < • • • < f3m{x) < Tm{x), Vx G [xo - S, Xo]. 
We can find ei = ei(S') such that for any e < ei and any e-generic x G [xq — S, xq] we have: 



34 BRANDON ROWEKAMP 

• UR^six) = ns{x) = m + l = ns{x^), and 

• the e-column of S^j.^.^ over xq consists of n{xQ) = m + 2 stacks. 
Theorem |4. 6| implies that 

'>^s,s{x) = ns{x) = nsix^) = m + 1 Vx E [xq — (5, xq — z^e] \ Ze. 

On the other hand, ns^ei^Q) = m + 2. Thus the interval [xq — ve, xq] must contain a jumping point 
of ns,e- 

Case 4. n(X([) < n(xo). This reduces the the previous case via the reflection 

9 ^ {-x,y) € M^. 

□ 



Remark 4.8. Theorem 4.6 states that the two functions n and coincide at points situated at a 



distance at least I'iS) pixels away from the jumping points of n. On the other hand, Theorem 4.7 
shows that, for a generic PL set, within ^{S) pixels from a jumping point of n there must be jumping 
points of □ 

Definition 4.9. Let 5 be a generic PL set in M?. 

(1) We will refer to the integer z^(5') as the noise range of S. 



(2) Let eQ{S) and £1(5") as defined in the Theorems |4 . 6 1 and 4.7 We set 

h{S) :=min(eo(5), ei{S)), 
and we will refer to it as the critical resolution of S. 

□ 



The next result explains the roles of the noise range and the critical resolution. 

Proposition 4.10. Let S be a generic PL set. If w > 2v{S), e < h{S) and [xq — we,xo + we] 
contains no jumping points ofn^, then n^{xo) = n(xo). 

Proof. Suppose that the interval [xq — we, xq + we] contains no jumping points of n^, yet ni;{xo) / 



Ti(xo). Then Theorem 4.6 implies that the interval 



w w 
[xo - —e,xo + —e\ 

contains a critical xi point of n. Then Theorem |4.7| implies that has a jumping point on the interval 

w w 
[xi - -e,xi + -e\ 

But since xi is at most -pixels from xq, this interval is contained within [xq — we,xo + we]. This 
contradicts the assumption that the interval contained no jumping points of and thus ne(xo) = 
n{xo). a 

Suppose that e < h{S), where 5 is a generic PL set. Then the theorems proven in this section up 
to this point imply that all the jumping points of are contained a fixed numbers of pixels from the 
jumping set of n. This simple observation, correctly implemented, will be the key to recovering the 
topology of S from the topology of its sufficiently fine pixelations. 

Consider the discontinuities of the function n^. They can only occur within i^iS) columns from 
a jumping point of n. We do not know what this integer is from the pixelation, but we know that it 
exists and it is independent of e. Therefore, we know that the noise range v{S) will eventually be less 
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than a properly chosen spread a{e) such that a{e) — )• oo as e — )• 0. Using the spread to estimate i^{S) 
will be a dramatic overestimation for small e. However since ea{e) — )• 0, if we declare any cycle 
which appears less than a{e) columns from a jumping point of ns{S) as a fake cycle, we will avoid 
declaring any real cycles as fake for small resolutions. 

The discontinuities ns{S) (for sufficiently small e) are obviously contained in the set 

{x G M : [x — ea{e),x + eo"(e)] contains a critical value of Ug}. 

We would like to consider this set to be the "noise portion" of S. This would mean that we could 
approximate S from PdS) over x-values outside of this region using the results from section 3. 
However, recall that in section 3 elementary sets were approximated by choosing upper and lower 
samples, which were chosen from the midpoints of pixels. This means that to use the methods from 
section 3 to approximate S we need our noise intervals to end at the middle of a pixel. With that in 
mind we define the noise interval: 

Definition 4.11. Let S be a generic PL set, and let 2:25 •• • , 2; at be the jumping points of ns,e- 
For each jumping point Xj let the noise interval Ijie) be the interval [aj, bj] where 

Xj E [aj,bj], 

aj{e),bj{e) £ -+eZ 

bj is the smallest number in ^ + such that (bj — xj) > 2ea{e) 

aj is the largest number in ^ + eZ such that (xj — bj) > 2ea{e) 
Define the set of noise intervals as 

= /i(e)U---U/jv(e). 

□ 

From the definition we see that A^ is a union of intervals. Furthermore for small e, each of these 
intervals will contain precisely one jumping point of ns- Thus A£(S') has as many components as 
the cardinality of Ss- In particular this implies that the measure of A^ is bounded from above by 
2|05|ecr(e). Since ea{e) vanishes as e — )• 0, this implies that A^ has vanishing measure as e — )• 0. 

Consider the set of rc-values which lie outside the noise. It is a finite union of intervals such that for 
all X which lie in these intervals, ni;{x) = n{x). The part of S situated above a each interval is either 
empty, or a disjoint union of regions of simply types. These types of regions can be approximated 
using the methods described in Section [3] Therefore, to complete the approximation of S, we need 
only describe how to deal with the noise intervals. 

The key observation is that the measure of Ag goes to zero as e — )■ with the properly chosen 
spread. This means that noise intervals make up a very small part of S, and so it will not be necessary 
to approximate them with as high of degree of accuracy as other parts of S. Indeed, we only seek to 
ensure that the noise intervals capture the correct homotopy type for small values of e. 

Since S is defined by a finite number of piecewise linear functions, for small enough e, the part of S 
above A^ consists of contractible connected components. Theorem |4. 1 [ implies that these components 
separate for small enough e. Therefore the easiest way to get the correct topology within the noise 
intervals is to simply cover each component of the pixelation above a noise interval with a rectangle, 
destroying any fake cycles (or holes) from Pe{S) (see Figure 11) 

Suppose that is a generic PL set with jump set 

3s = {xi < ■ ■ ■ < xn}- 
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Figure 1 1 . The approximation of a noise interval. The red rectangles indicate the 
rectangles that form the approximation. Since the intersection of a noise interval and 
a regular interval always lie on the centers of a column, the edges of the rectangles 
will intersect approximations from the regular intervals. 



We describe below an the algorithm to create a PL approximation SgV for S using only information 
from Pe{S) is as follows (this algorithm is restated in a more concrete way in Appendix [P]). 

Algorithm 4.12. (1) Choose a spread a such that ea{e)'^ — )• oo and ea{e) — )• as e — )• 0. 

(2) Let Ae be the noise intervals of P£{S). 

(3) Define Ol^ to be the closure of M \ A^. We call the regular set, and its intervals are called 
regular intervals. 

(4) For the bounded regular interval / € do the following: 

(a) For each connected component of Pe{S) H (/ x M), the part of P£{S) over the regu- 
lar regular intervals / choose compatible upper and lower profiles 11+ and 11^ on that 
component. 

(b) For each connected component, choose compatible upper and lower samples Ef and S~ 
with spread a. 

(c) For each connected component, take the PL approximation of the upper and lower sam- 
ples. 

(d) Add each PL approximation to the approximation of S. 

(5) For each regular point x denote by Se{x) the part of Se over x, i.e., 

Se{x) :=Sen{x} xM. 

(6) For a noise interval Ij{e) we observe first that Corollary |4.2| implies that the number of con- 
nected components of PeiS) over Ij (e) is equal to the number of components of S over Ij (e). 
Denote by Qj{e) the set of connected components of Pe{S) over Ij{e) 

(a) For every C £ Qjie) of Pe(5) we set 

C/c :=max{y; {x,y)GC}, Lc := min{y; (x,y)eC}. 

(b) Let P^ denote the pixel whose top boundary lies on the line {y = Uc} and intersects the 
line {x = Ss{aj{e))}. Similarly let Pjj denote the pixel whose top boundary lies on the 
line {y = Uc} and intersects the line {x = let P£ denote the pixel whose 
bottom boundary lies on the line {y = Lc} and intersects the line {x = 5e(aj(e))} and 
P^ denote the pixel whose bottom boundary lies on the line {y = Lc} and intersects the 
line {x = Se{bj{e))}. (These pixels lie on the comers of a rectangle which bounds the 



majority of C, see Figure 11) 
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(c) Let Uq indicate the center of and similarly let Uq, Uf^ and indicate the centers of 

Pfj, PI and Pf . 

(d) Denote by 5*0 (e) the convex quadrilateral with vertices 

BJ{e) = {h,{e),L'c), := f/^), 

A-{e) = {a,{e),L-c), A+{e) := {aj{e),U^) 
(7) To the set Se constructed at (4) add the noise union 

N 

^-=U U ^cie). 
j=iceejie) 

This final set will be piecewise Unear by construction, and will be a good approximation of the 
set. □ 
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5. The Main Result 

Consider a generic PL subset S of the Euclidean plane E = M^. In section|4]we created Algorithm 



4.12 which creates a PL approximation 5^ of S. We wish to state that this approximation converges 
in a good way to S. An appropriate language to state this condition is that of normal cycles. 

For any compact PL subset X C E we denote by A^"**- its normal cycle, EUTl. For the reader's 
convenience we have included in Appendix [C] a brief survey of the basic properties of the normal 
cycle. In particular Appendix [C] demonstrates how to recover both the Euler characteristic and the 
perimeter of a set using only calculations on the normal cycle. Using similar techniques we can extract 
other geometric and topological information about a set using only its normal cycle. Therefore weak 
convergence in normal cycles implies convergence of a great deal of important information, and is an 
appropriate condition for "good convergence." 



In the final theorem we see that the approximation Ss created by Algorithm 4.12 does satisfy this 
type of convergence. 



Theorem 5.1. Suppose S is a generic PL subset ofK, a{e) is a sprea d function satisfying {3.1a ) and 



Ss is the PL approximation of S constructed via the Algorithm 4.12 Then the normal cycle of 
Ss converges weakly to the normal cycle of S. 

Proof. Let us first outline the strategy. Recall that V5 indicates the set of vertices of S and denote 
by its projection on the x-axis. Since S is generic the projection E 9 (x, y) 1— > x G R induces a 
bijection V5 — )• X5. The jump set is contained in X5. We will refer to the vertices that project in 
05 as essential vertices. The other vertices will be called inessential. 

For every c G cJs and any e > small we denote by Is{c) the noise interval that contains c. We set 
/o(c) := {c}. Fix eo > such that, 

dist(A£(c), Ae(c)) < ^dist(c,c'), \lc,c^ds, < e < eo- 
For c G and any < e < eo we denote by !N"e(c) the noise strip 

?^,(c) := {(x,y) gM2; x G 4(c) }. 
Note that No(c) is the vertical line {x = c}. Finally we set 

X, = IJ N,(c), 'Rs ■■= \ N„ VO < e < eo- 

For uniformity, we set Sq := S. 

Observe that there exists < ei < eo such that, if e G [0, ei], the following hold. 

• Any component of n Sq is contained in a unique component of 3sfe n S'g. 

• If C is a connected component of 3?^ n Ss, then the closure of C intersects exactly two 
connected components of Jsf^ n S^. 

For any e G [0, ei] we construct a graph as follows. The vertex set of consists of the 
connected components of !Ne n S'g. The edges are the connected components of 3?^ n We have a 
well defined map 

: Vn ^ Vs 



that associates to the a component C of ?\fo n 5 the unique component of Ss fl l^s that contains C. 
which This is easily seen to be a bijection. Moreover, it induces an isomorphism of graphs, i.e., the 
vertices v, v' are adjacent in Fq if and only if the vertices ^'e(i') and ^'^(ij') are adjacent in Fg. For 
any t; G Vq we denote by the set of edges of Fq that are adjacent to v. 
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Remark 5.2. The graph Fq is known in the literature as the Reeb graph of the (stratified) Morse 
function 

For more information about this concept we refer to pi, VI.4], or the original source ll22l . □ 

For any vertex G "Vq and e G [0, ei] we denote by C„ ^ the component of n Se corresponding 
to '^^{v). Similarly, for any edge e = [v, v'] of Fq and any e G [0, ei] we denote by Ce,e the closure 
of the component of 3?^ n corresponding to the edge \^e{'i'), '^e{v')]- 

Lemma 5.3. For any e G [0, ei] we have 

JSfS, ^ J^C^,, _^ ^ jyCe,. - ^ ^ ]^C^,er^C^,e _ (5.1) 

iieVo eeeo v&Vo sgE-u 

Proof. The proof use the inclusion-exclusion principle, i.e., the equality 

for any compact PL subsets X,Y C M?. 
Note that we have a decomposition 

Se = I \J C^,e] 'J I [j Ce,e ) • (5.2) 

We need to discuss separately the cases e > and e = 0. 
1. Assume that e G (0, ei]. In this case we have 

C„,. n C„/,, = = Ce,e n C7e.,„ \/v^v', e^e'. (5.3) 



The equality (5.1 1 now follows from inclusion-exclusion principle applied to the decomposition (5.2 1 
satisfying the overlap conditions (5.3 1. 

2. Assume that e = 0. In this case the overlap conditions are more complicated. We have 

Cvfi n (7„/,o = 0, Vtj / v', (5.4a) 
Ce,o n Ce/,0 = ^ e n e' = 0, (5.4b) 
where the condition e n e' = signifies that the edges e and e' have no vertex in common. Moreover, 

n Cefi = C,,o, G Vo, Ac E^. (5.5) 

eeA 

Using (5.2 1, ( |5.4a| ), ( |5.4b i, < \5.5} and the inclusion-exclusion principle we deduce 

i>eVo eG£o rieVo eSS^ 

+ ^ ^ (_l)l^l+iArf^.n(agA)Ce,o + ^ ^ (_l)l^lArneeAC«,o 

_ 'y ^ _j_ ^ ^ J^l^efi _ 'y ^ ^ ^ jyC'ij.onCe.o 

weVo ee£o vGVo ei^E^ 

+ E I E ( (-1)1^1+^ + (-1)1^1 )) 

^ V ' 

=0 

□ 
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Theorem 5. 1 is now an immediate consequence of the following result. 



Lemma 5.4. For any v G Vq, e ^ Zqwc have 

limTV^'''^ = TV^-'O, (5.6a) 

e\0 

lim A^'^'='- = 7V^<='0, (5.6b) 
lim ArC->-nc.,, ^ ^Ce,onc„,o^ ^g ^^,) 

e\0 

where the limits are understood in the sense of weak topology on the space of currents. 

The proof of this lemma relies on the General Convergence Theorem proved by Joseph Fu in ||6l. 

Theorem 5.5 (Approximation Theorem). Suppose S is a PL subset of the plane and for each e G 
[0, e\\ Ss is a PL subset of the plane with the following properties. 

(1) There is a compact set K cM? which contains each Se. 

(2) There is a M ^M. such that 

mass(A^'^=) <M, Ve. 

(3) For almost every ^ G Hom(M^, M) and almost every c G M we have 

i\mx{Sen{c>c}) = x{Sn{c>c}) 

Then N^^ converges to N''^ as e ^ weakly and in the flat norm. □ 



Proving ( 5.6a i-( 5.6c i will involve proving each of the three conditions in the Approximation The- 
orem. We begin with an easy consequence of the Approximation Theorem that will be very useful 
in the sequel. First, let us define a convex polygon in the plane to be the convex hull of a finite set. 
Note that this definition allows for degenerate polygons such as line segments or points. We define 
the perimeter of a segment to be twice its length. For 2-dimensional polygons the perimeter is defined 
in the usual way. We will denote by L{P) the perimeter of a polygon. 

Lemma 5.6. Suppose {Se)e>o is a family of convex polygons in the plane that converge in the Haus- 
dorff metric to a convex polygon S. Then N^"" converges weakly to as e ^ 0. 

Proof. We argue by proving the conditions of Fu's Theorem. Observe first that there exists R > 
such that 

dist(5£,5) < R, Ve 

and thus the condition (1) of the Approximation Theorem. The computations of [20, Chap. 23] show 
that mass of the normal cycle of a convex polygon P is equal to 2ir + L{P). From Hadwiger's 
characterization theorem lITSl Thm. 9.1.1] we deduce that 

limL(5e) = Lis) 

and thus condition (2) is also satisfied. 

Therefore we must show that for almost every ^ G Hom(M^, M) and almost every c G M we have 

iimx(5enU<c}) = x(5nU<c}) 

£\0 

Note that S and each 5^ are all convex subsets of the plane. Therefore any intersection with a half- 
plane will either be empty or be a contractible set. Therefore to prove the convergence of Euler 
characteristic on half-planes we need only prove that a half plane H will only intersect Se for small e 
if and only if it intersects S. This is true since H n Se converges in the Hausdorff metric to H n S. □ 
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In several places in the remainder of the proof of Lemma 5.4 it will be convenient to discuss the 



maximum rate at which S can increase or decrease. With that in mind we set: 

a = sup |/3| 

/3 is a slope of an edge of S 



Proof of ( 5.6a| ). Note that each C„^o is a connected component of a subset of a line. Therefore 



each Cyfi is either a point or a line segment. By construction of the approximation Sg, each C„_£ 
is a rectangle which contains C„^o- Therefore to show that C„ ^ converges to C^fl in the Hausdorff 
metric we need only show that its width vanishes and its height converges to the height of C„^o- By 
construction the width of C^^e is 0{ea{e)), and since cr is a spread this width vanishes as e — )• 0. 



Furthermore Proposition 2.3 implies that the difference in the height of C^^e and C^fl is 0{ae + 
aea{e)), where both as and aea{e) vanish as e — )• 0. Therefore C^,_e converge to C^,_o in the 
Hausdorff metric, and Lemma [5^ implies that 



lim iV^"-' = iV'^^'O 

e\0 



Proof of ( 5.6c| ). The set C„ ^ n Ce,e is a line segment on the edge of the rectangle Cv,e- We have 



already proven that C^^e converges to C„ o in the Hausdorff metric. Since 

Cvfl n Ce,0 = Cvfl 

to prove that Cv,£ H Ce,e converges to Cv,o n Ce,o in the Hausdorff metric we need only prove that 

lim sup inf dist(x,y) = 0. 

That is to say, we must show that the maximum distance from a point in C„ o to the set Ce,s H 
becomes arbitrarily small as e — )• 0. But note that Ce,e H Cv,e is an interval which lies on the edge 
of a noise interval. Its distance to a point of Cvfl in j;-coordinates is at most the width of the noise 
interval. Its distance to a point of Cvfl in y-coordinates is equal to the change of slope of S over the 
noise interval, together with the error from S to the approximation Se. Thus (once again using the 



e 



fact that a noise interval has width proportional to ea{e) together with Proposition 2.3 1 we see that 



sup inf dist(x, y) = 0{ea{e) + ae) 

which implies that the Ce,s H Cv,e converges to Cvfl = Cefl H Cvfl in the Hausdorff metric. So 
Lemma [53] implies that 

lim jV'^'^'''""^'"'^ = ]\[^'''0<^^v,o 
e\0 



Proof of (5.6b ;) To prove this hmit we make further use the inclusion-exclusion principle. 

Fix an edge e. Let the set Xg o indicate the projection of vertices of Cefl to the x-axis (i.e. the 
subset of Xs which contains only projections of vertices in Cefl)- Write 

•^e,0 = { Xq, Xi, . . . , Xn } , 

where xq, xi, . . . , x„ are arranged in increasing order. For each integer 1 < i < n we define the set 

Vx.fl ■■= Cefl n{x = Xi} 

and for each integer 1 < i < re we define the set 

Rx.fl ■■= Cefl n {Xi^i < X < Xi}. 

Note that Rx,fl n Rx,+,fl = 14„o- 
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Let Xe e be the projection of the vertices of Ce,e to the x-axis. Note that by construction of Se each 
point in Xe,e will have exactly two vertices of Ce,e map onto it. For each e write 

■^e,e = { ^0,e^ • • • > ^m{e),£ } 

arranged in increasing order and where m(e) is an integer depending on e. 

Where jj ^ is defined to be the largest integer such that Xj^ < Xi and conversely ki^e is defined to 
the smallest integer such that x^- > Xj. The set ^ can be thought of as the part of Ce,£ between 
the line segments which cross over Xi (see Figure 12 1. Note that if the critical value Xi belongs to 
Xe,e, then the set T^. ^ will be a hexagon, otherwise it will be a quadrilateral. 



Figure 12. In this picture ttie grey region is Ce,e, the red line is {x = Xi} and the 
blue region is Vxi^e- 



Similarly, for each e > and each integer \ <i <n define the set 

where jj e and ki^^ are defined as before. Then the set Rxi,£ is the part of Ce,e which lies between 
and K„e- 

Dividing Ce,o among the sets Rx^,o and recalling that these sets overlap on Vx^fl we see that 

n n—1 

AtC-.o = j^R.,A-> _ j^v^.fi (5 7) 

i=l 1=1 

Dividing Ce,e among the sets Rxi.e and the sets ^ we see that 



n-l 



N 



n—1 n—1 



(5.8) 



i=l 



1=1 



i=l 



Although (5.8l looks considerably more complicated than (5.7| note that (5.8l will converge to (5.7i 
if the following three equations are true: 



£\0 



(5.9a) 
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Figure 13. The region T4-,e fiay be a hexagon if the line {x = Xi} occurs at a 



vertex of C, 



e,e- 



lim AT-^^i.^'^^^*.'^ = lim AT^^i'^*^-^^*- 
We will prove these equations in reverse order. 



(5.9b) 
(5.9c) 



Proof of ( 5.9c| ). We start by proving 



lim N^-^■^^^-. 
e\0 



Note that the intersection of Rxi,e and Vx^^e is a line segment. The set T^. o is also a line segment, 



so to prove the convergence of normal cycles we can use Lemma 5.6 For sufficiently small e we 



can assume that no vertex of S occurs in between Rx,,e H 14- ,e and Vx^fl. Then the top vertex of 
R'Xi,e n T4i,£ will differ in x-coordinate from the top vertex of I4.,o by at most the width of a line 
segment in Ce,e and differ in y-coordinate by a number proportional to the width of the line segment. 
This means that the distance from the top vertex of Rxi,e H 14, ,e to the top vertex of T4. o is 0{ea{e)). 
The same is true of the bottom vertices. Therefore Lemma [53] implies that 

The above reasoning can be repeated to show that 



lim N^^^'^^^'^i^ 

e\0 



N 



which completes the proof of (5.9c i. 

Proof of ( 5.9b[ ). We once again make use of the inclusion-exclusion principle satisfied by the Normal 
Cycle. Divide the set 14- ^ into the following two sets: 



Then we have 



((00,Xi 



1/0 



{[Xi,00) X R)nVx^,e 



SO that (5.9b I will hold if iV^--s A^"^--^ and N 



all converge to N^^'i-^. 



Note that the set V^. ^ n Vy. ^ is a line segment. If the top and bottom vertex of 



sample points of the approximation, then Proposition 2.3 implies that the Hausdorff distance from 



x„£ n V^^, are 



14. ^ n 14^ £ to 14., is 0{a£). If the top and bottom vertices of the line segment are not sample 
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points of the approximation, note that they differ in width from a sample point by at most ea{e). 
Therefore in all cases the Hausdorff distance from V^-^i; H V^.^^ to Vx^fl is 0{ae + aea{e)). This 
vanishes as e — )• 0, so that Lemma [53] implies that 



£\0 



Note that each set V^^. ^ and V^-^^ must be a convex quadrilateral. The width of their union Vx^^e is 
0{ea{e)) by construction. This implies that the widths of each of 1^9. ^ and Vx. ^ vanish as e — )• 0. 
Therefore each set converges to 1^9 ^ n Vx. ^ in Hausdorff distance as e — )• 0. This implies 



I/O 

lim N 

e\0 



T/l 

lim N 

e\0 



completing the proof of (5.9b i. 
Proof of ( 5.9a) ). The major difficulty in this case is that the regions Rxi,£ need not be convex polygons 



so we cannot apply Lemma 5.6 Thus we must return to proving the three conditions of Theorem 5.5 
directly. 

Note that each Rx^.e can be considered the PL approximation with spread (j{e) of the region Rxi,o- 



The region Rxi,o is simply an elementary set which lies between two line segments. Theorem 2.11 
implies that the distance of the boundary of Rxi,£ to the boundary of Rxi,o is 0{e+ea{e) + (e(T(e))^). 
Therefore the sets Rx,,e must all be contained in some compact subset of M?. 



Since Rxi,£ is the PL approximation of an elementary set with spread a{e), Corollary 3.7 implies 



that the total curvature of its boundary converges to the total curvature of the boundary of Rxi,o- The 
computations of ll20l Chap. 23] show that the mass of a normal cycle of a polygon is equal to the 
perimeter of the polygon plus the total curvature of its boundary. Therefore the mass of A^^^i -^ cannot 
increase to infinity. 



The final condition to prove from Theorem 5.5 is that the Euler characteristic of intersection of 
Rx,,£ with any generic half -plane converges to the Euler characteristic of the intersection of Rx,,o 
with the same half-plane. To prove this we will need a technical lemma whose proof we will defer 
until later. 




Figure 14. The situation of Lemma 5. 7 The red line only intersects the boundary 
of the blue approximation twice. 
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Lemma 5.7. Let S be a single non-vertical line segment and for each e > let an upper or lower 
sample of S with spread a{e) and for every e let fs indicate the PL-approximation of S determined 
by Ee. 

Let ^ G Hom(M2, M) and suppose that the line = c} is not parallel to S and intersects S. Then 
there exists a e* > depending only on S and ^ such that for < e < Eq, the line = c} intersects 
the graph of exactly once. □ 

Given this lemma, fix a c G M and a ^ G Hom(M^, M) such that the line = c} is not parallel to 
a line segment of Rx^,o- We distinguish two cases. 



A. The line = c} does not intersect Rxi,o- Theorem 2.11 implies that for all sufficiently small 
e, the line = c} will not intersect Rxi,s- Since Rxi,e converges to Rxi,o in Hausdorff measure, 
the half -plane will contain Rx^,e for small e if and only if it contains Rxi,o- The sets Rx,,o and Rxi,e 
both must be contractible since they are regions which lie between the graphs of two piecewise linear 
functions which do not intersect. Therefore in this case 

lim > c} n Rx,,e) = Xi{^ > c} n Rx,,o). 

B. The line = c} does intersect the set Rxi,o- Since Rxi,o is convex, we clearly have 

x({e > c} n Rx,,o) = 1. 



Let e G (0, e*{Rxifi)), where e*{Rx^fi) is the constant guaranteed by Lemma 5.7 Then the desired 



convergence in Euler characteristic is an immediate consequence of the following lemma. 

Lemma 5.8. (a) For any d > c the intersection between the line = d} and Rxi,e either empty 
or a closed segment (possibly degenerate). 

(b)x[{^>c}r^Rx,,e) = l. 



Figure 1 5 . For a line to intersect Rx^ ,e in a non-connected fashion it must intersect 
either the upper or the lower boundary multiple times. 



Proof, (a) Assume that the line = d} is the graph of the linear function l{x) = mx + b. The 
boundary of Rxi,e has four component: a top, a bottom and two vertical side components. The top 
and boundary components are graphs of PL functions Ts{x) and respectively Be{x) (see Figure 15 1. 
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Lemma 5.7 implies that each line = d} intersects each of these components at most once so the 
intersection between the line = d} and the boundary of Rxi,e consists of at most four points. 

We consider the set = d} Ci Rx,,£ which is the union of some number of line segments. Let 
pi, P2, ■■■ Pn, n < 'i,be, the endpoints of these line segments arranged increasingly according to their 
x-coordinates. We claim that n < 2 so that the intersection between = d} and Rxi,s is a (possibly 
degenerate) line segment. To prove this we argue by contradiction. 

Suppose that the intersection consists of at least three points pi,P2jP3- Let {tj,Uj) be the coordi- 
nates of Pj, I < j < 3, so that ti < t2 < ts- Think of the line = d} as the trajectory of a particle 
moving in the plane with constant velocity 

X = t, y = £{t) = mt + b. 

The particle first enters Rxi,e at the moment t = ti and exits at the moment t = t2- Note that the 
corresponding exit point p2 cannot be on either one of the vertical components of the boundary of 
Rx^,£ so it must be either on the top part or on the bottom of the boundary. Assume that p2 is on the 
top part of the boundary. (The case when p2 is on the bottom is dealt with in a similar fashion.) 

Lemma [5^ implies that the particle will never intersect the top part of the boundary ever again so 
that 

i{t) > Te{t) > Beit), Vt > t2. 

This implies that at the moment ts when the particle enters Rxi,e again, it must do so through one of 
the vertical portions of the boundary of Rxi,e- On the other hand at that moment we have 

U3=eih)>Teit3)>Beih). 

so the point p^ cannot lie on the vertical segment 

{{t3,y);B,{t3)<y<T,{t3)} 

which the rightmost part of the boundary of Rxi,e- 
(b) Set for simplicity 

Rl^,l:= {C>c}nR,,,s. 
Observe that the set Ic := ^ (-Sl^e ) C M is connected because ^(Rxi^e) is a closed interval and 

Ic = C{Rx,,e)n[c,00). 

The map ^ defines a continuous surjection ^ : -R|~e — )• Ic with contractible fibers. The Vietoris-Biegle 
theorem, [23^ Thm.l5, Chap. 6, Sec.9] implies that i?|~e has the same coh □ 



To complete the proof of ( |5.9a i (and consequently also prove (5.6bi, Lemma 5.4 and the overall 



Theoren jSTT] ) we need only prove Lemma 5.7 



Proof of Lemma 5.7 Let a be the slope of 5. Suppose that ^(xi, X2) = ^1X1 + ^2X2- Then 

since S is not parallel to = c}. The function is differentiable every where except at a finite 
number of points. Theorem |2. 11 implies that there exists an eq > such that for < e < Eq, 

is either strictly positive or strictly negative for all x in the domain of fe since /^(x) must be close to 
aand,^i+6a 7^ 0. Weset^|/^(x) = Cix + 6/(x). The notation indicates that ^|/^ is the restriction 
of ^ to the graph of /g. Note that ^|/^ is Lipschitz and 

^(x) = 6 + 6/;(x), 
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for all but finitely many x. 

Now suppose that {i = c\ intersects the graph of at points (xq, yo) and (xi, yi). Then note ^| 
is equal to c for xq and xi. Thus we have, [24, Prop. 11.12] 



dx 



dx 



{Ci + ^2fUx))dx. 



XO 



But ^1 + ^2/e(2;) is either strictly positive or strictly negative. Therefore the integral 

del/. 



XO 



dx 



-dx 



cannot be 0, contradicting our assumption that = c} intersects the graph of fe at two points. This 
completes the proof of Lemma [5^ □ 
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Appendix A. The Farey Series and Holes in Pixelations 

An unintuitive feature of pixelations is that tiiey do not preserve iiomotopy type, even for small 
values of e. The Figure 10 in Section |4] gives an example of a contractible set whose e-pixelations 
always contain a cycle. In this appendix we wish to expand on the nature of these "holes" as well as 
show that an arbitrary number of such holes can appear. 

We start by defining a convenient class of sets. 

Definition A.l. The angle with slopes (a, where /3 > a is the set 

A{^,a) := {{x,y) eR^ : x>0, (y - ax)(y - /3x) = }. 



Angles are convenient to work with in the context of pixelations since they have nice self similarity 
properties. Note that 

PeiA{l3,a))=ePi{A{l3,a)) 
where the last equality holds because A{P, a) does not change under rescalings centered at the origin. 
This means that different e-pixelations are simply retractions or expansions of any other e-pixelation, 
and are thus topologically equivalent. Therefore if a cycle appears in any e-pixelation of an angle, 
it will appear in every e-pixelation of an angle. We refer to these false cycles in the e-pixelations as 
holes. 

Figure 16 is an example of a pixelation of an angle that has a hole, and in particular, the pixelation 
is not contractible, while an angle plainly is. 



Figure 16. The pixelation of the angle ^4(1,^) contains 1 hole in the 7th column 
from the y-axis. 

When varying the slopes of the angle A{I3, a) several things become apparent. First, for a hole 
to appear either both slopes must be positive or both slopes must be negative. Second, the two lines 
must have fairly close slopes for holes to appear. Finally, more holes will tend to appear when one of 
the slopes is close to 1. 

It is beyond the scope of this paper to classify the behavior of holes for general angles. However, 
an easier situation occurs when the two slopes of the angle are adjacent members of a Farey series. 

The n-th Farey series F„ is the increasing finite sequence consisting of the rational numbers be- 
tween and 1 which have denominator of n or less when written in lowest terms. For example, 

0111231 
^~ 1'4'3'2'3'4' 1 
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An important property of the Farey sequences is that if | and ^ are consecutive terms in a Farey 
series, then be — ad = 1, |[T2l Chap.III]. Using this property, we can prove the following fact. 

Proposition A.2. Suppose ^ and ^ are consecutive terms in a Farey series. Denote by g{a, b; c, d) 
the number of holes of the e-pixelation 

and by Xa^b,c4 of lattice points {p, q) G such that 

0<p<d, 
a < q < min(2a, b), 
bd + bc + b + d<l+pb + qd<2bd + bc + b + d. 

Then 

g{a, b; c, d) := ad - #Xafi,c,d 



Proof. As noted in the earlier discussion, any e-pixelation of the angle is topologically equivalent to 
Pi (j4 ^)). Therefore we need only to consider the case e = 1. We set 

c a 
mi := mo := -. 
d b 

For j = 0, 1 we denote by ij the half-line 

y = rrijX, x > 0. 
For each i G Z>o we define (using the notations in Definition |2.1[ ) 

'[mi(i-l)J, mi(i-l)0Z 



Li 



2i-l 



2i - 1 



mi(i — 1) — 1, mi(i — 1) G 



(A.2) 



(A.3) 



niQi + 1, niQi G Z 

In other words, Ui indicates the lowest y-value of the 1-pixelation of the half-hne £i within the i-th 
column [i — x M, and Lj indicates the highest y-value of the 1-pixelation of the half-line £q in 
the same column. Set Di := Ui — Li and observe that ( |A.2 ) and ( A.3| ) imply 



Di < (mi — mo)« 



bd' 



Note also that 



More generally, ( A.4| ) implies that 



D, = -2, D2£{-l,-2}. 



Di < 0, yi = l,...,bd. 



(A.4) 
(A.5) 
(A.6) 



The quantities Ui and Li are useful because they allow us to easily state when there is a gap in a 
column or on the edge of a column. The i-th columns of the pixelations of the upper and lower 
half-lines overlap if and only if 

Ui < L, 
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and the two pixelations overlap on the right side of the i-th column if and only if 

Ui < Lj+i 

A hole is a gap which is closed on both sides. In particular, it is closed on the right side. This means 
in each hole there is a column where the pixelations of £q and ii do not overlap on the interior of the 
column, but do overlap on the right side of the column. In terms of Ui and Li this is equivalent to the 
condition 

Li<U< Li+i (A.7) 
Note that this condition only occurs when Lj+i > Lj. Since ^ < 1 we have 

Li+i < Li + 1 

This means that the condition L^+i > Lj is equivalent to the condition L^+i = Lj + 1. 

Suppose that Ui — Li = 2. Then it is impossible that Ui < Lj_|_i since that would imply that the 
lower line increased by at least the length of a pixel in the amount of time that the upper line increased 
by less than the length of a pixel. Therefore for (A.7 1 to occur, it must be true that Ui — Li = 1. 
Combining this fact with the above discussion about Li and Lj+i, we find that (A.7 1 is equivalent to 
the condition 

Ui -Li = l and Li+i = Li + l. (A.8) 
Note that Li increases by c every d columns since ii has slope ^ . Hence 

Ui+d = Ui + c, 

and similarly 

Li+b = Li + a. 

Therefore 

Ui+bd - Li+bd = Ui + bc- {Li + ad) = U - Li + {be - ad) = C/j - Lj + 1, 

so that 

D^+bd = A + 1. (A.9) 

We set 

ki := min{A; E Z>o; = 1 }• 
We know that ki exists since (A.9 1 and ( A.5 1 imply that Di takes on every integer value greater than 
or equal to 1. Note that ( A.4 1 implies that, if ki exists, then ki > bd. To finish the proof we will need 
a variety of facts about ki. 

Lemma A.3. (a) ki := {a + l)d + b + 1 + bd = bd + be + b + d. 

(b) For any i £ [ki,ki + bd) CiZwe have 

< A < 2 (A.IO) 

(c) For any i > ki + bdwe have Di > 0. 



Proof, (a) The inequality ( A.4 1 implies that ki > bd so that 

ko := ki-bd> 

and 

Dk, = 0. 

Moreover, the minimality of ki coupled with ( |A.9 1 implies that 

ko = min{i > 0; Di = 0}. 
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If = then there exists a positive integer £ such that 

£ - I < mo{i - 1) < ruQi < i. (A.lla) 

£ <mi{i-l) <£ + l (A.llb) 

When these conditions are satisfied we have Ui = Li = £. The lattice point (i — 1, ^) is in the interior 
of the angle A(mi, mo) spanned by the vectors 

ui = {b, a), U2 = {d, c). 

Since be — ad = 1 we deduce that the vectors ui, U2 form an integral basis of the lattice Z^. We 
deduce that there exist two positive integers p, q such that 

{i — !,£)= pui + qu2, i.e., i — 1 = pb + qd, £ = pa + qc. 

Observe that 

/ . , , a , , ,s ad be — I Q „ Q 

niQii - I) = -{pb + qd) = ap + q— = ap + q — ; — =ap + qe --=£--. 

b b b b b 

Hence ( A. lla| ) is satisfied if and only if 

a<q<b. (A.12) 

On the other hand, 

c be ad 1 P P 

mi{i - 1) = -Apb + qd) = p— + qe = p \- qe = pa + qe + - = £ + -. 

d d d d d 

Hence ( |A. llb| ) is satisfied when 

Q<p<d. (A.13) 

From the equality 

i = 1 + pb + qd 

we deduce that 

min{i > 0; A = 0} = min{l + pb + qd; < p < d, a<q<b} = l + b+{a+ l)d. 
This proves (a). 

(b) To prove the upper estimate in ( A. 10| ) we use ( |A.4 1. For i < ki + bd v/e have 

ki (a + l)d + b+l a + l 6+1 

''•<'^^ = '^- M = l + ^ + l^<^- 

Thus Di < 2 for all i G [ki,ki + bd) n Z. \f Di = 2 for some i in this range then Di^M = 1, 
contradicting the minimality of ki . 

To prove the lower estimate part of (A. 10 1 we recall that 

ko = ki — bd = min{ k > 0; = }. 

Let us observe that 

Di > -1, Mi > ko. 

Indeed the inequality Di < —1 takes place only if niQi > ^ G Z and mi(i — 1) < £. This implies 

% c 

niiii — 1) < rriQi <;=^ ^ < -; <^=^ i < be = ad + 1 < k^. 
bd d 



Hence, for any i > kQ + bd v/e have Di > proving part (b). Part (c) follows from (b) and ( A.9 1. □ 
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This lemma narrows down the possible locations of holes. In particular, property (c) implies that 
holes can only occur in columns in the interval [A;i , ki + bd). Recall that for a hole to appear, we need 
Lj+i = Li + I and Di = I. In most cases when Lj increases in this interval Di will also be equal to 
1. However, if {7j+i = C/j + 1, Lj+i = Lj + 1 and Di = 0, then no hole appears. We want to count 
these instances where an increase in Li does not indicate the presence of a hole through the use of the 
set Xa,b,c,d- We do this through a series of lemmas. 

Define a map 

: Xa,b,c,d ^, {p, q)^l+pb + qd. 
Lemma A.4. The map <I> is injective. 

Proof. Indeed if <I>(p, q) = ^{p' , q') then {p — p')h= {q' — q)d. Since b and d are coprime we deduce 
that [p — p') must be a multiple of d. Since < p,p' < d this happens if and only if p = p'. It 
automatically follows that q = q' . □ 

We denote by Ia,b,c,d the range of <I>. From the definition of Xa^b,c4 we deduce that 

Ia,b,c4 C [ki,ki + bd). 

Lemma A.5. The following statements are equivalent. 

(a) i G [/ci, ki + bd) OZ, Ui = Li and Lj+i = Lj. 

(b) i £ Ia,b,c,d 

Proof. Let i £ [ki,ki + bd) Pi Z such that Lj+i = Lj + 1 and Di = 0. Using the notation and the 
terminology employed in the proof of Lemma |A.3|^a) we deduce that the condition Di = holds if 



and only if that there exist positive an integer i such that ( A. 1 1 a i and ( A. 1 lb i hold. Moreover 

i = 1 + pb + qd., 



where q and p are constrained by (A.12i and (A.13l. Once the condition Dj = is satisfied the 
condition L^+i = L j + 1 = + 1 is equivalent to 

I < mo{i + !)<£+!. 
moil + 1) = mo{i - 1) + 2mo =£ + ^ > i. 



Hence a <q< min(6, 2a), < p < d, and i = ^{p, q), {p, q) G Xa,b,c,d, i.e., i G Ia,b,c,d- 



□ 



Lemma A.3 c) shows that the conditions of (A.8l cannot be met in the {k + bd)-th column and 
beyond. Therefore, the right edge of any hole must occur between the A;-th column and the {k + bd)- 
th column. Therefore the pixelation of the angle contains a hole for each time that Ui — Li = 1 and 
Lj increases between the k-th and {k + bd)-th column (including k but not {k + bd).) and This is 
a total of bd columns, and Lj increases by a every b columns. Hence Lj will increase ad times in 
this range. On the other hand if z G Ia,b,c,d the Lj increases yet Ui = Li. Therefore (^A ^)) 
contains ad — #Ia,b,c,d holes. The conclusion of the proposition follows from the injectivity of <I>. □ 



Remark A.6. The proof of Proposition A.2 gives us a simple algorithm for computing g{a, b; c, d) 



that is easily implementable numerically. More precisely we have 



(o, 6; c, d) 



ki+bd-l 
i=ki 



Li){Li^ 



Li 



□ 
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Corollary A.7. 



g{n,2n + 1; 1,2) = 2n. 



□ 



Proof. In this case the set X„_2n+i,i,2 is empty and the above equality follows immediately from 
Proposition A.2 □ 



The above corollary shows that the pixelation of even a very simple set can have a complicated 
homotopy type. 



Example A.8. Let us consider the following situation 

c 1 a 

Then 



i+2 



■j+5 



Lj + 2, 



A; 



-10 



A; + 1. 



Using (A.2 1, ( A.3 1 and the above equalities we get 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Ui 


-1 








1 


1 


2 


2 


3 


3 


4 


4 


Li 


1 


1 


2 


2 


3 


3 


4 


4 


4 


5 


5 




-2 


-1 


-2 


-1 


-2 


-1 


-2 


-1 


-1 


-1 


-1 



Using this table and ( A.9 1 we deduce that k 

Ui+20 

we obtain the following table 



22. Next, using the above table and the equalities 

Ui + 10, Li+20 = Li + 8 



i 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


Ui 


10 


10 


11 


11 


12 


12 


13 


13 


14 


14 


Li 


9 


10 


10 


11 


11 


12 


12 


12 


13 


13 


Di 


1 





1 





1 


1 


1 


1 


1 


1 



In this case there are five holes. Two holes are immediately apparent, since in the 23rd and 25th 
column a distance of 1 is followed by a distance of 0. However note that in the 26th column we have 
£^26 of 1 and then C/26 = ^27- This indicates that a hole is formed merely because the pixelations of 



the upper and lower lines touch at a comer (compare to the case of Figure 16 ) Similar holes exist in 
the 27th column and the 29th column. No holes exist past this point, since Lemma A.3 implies that 
past this point Di > 0. □ 



Appendix B. Subanalytic Currents 

In this appendix we gather without proofs a few facts about the subanalytic currents introduced by 
R. Hardt in ifTOl [Til . Our terminology concerning currents closely follows that of Federer [5 ] (see 
also the more accessible |fT6','20|). However, we changed some notations to better resemble notations 
used in algebraic topology. First we need to define the subanalytic sets. 

An M.-structure is a collection § = { S" }^^-^, S" C T(M"), with the following properties. 

El.: S" contains all the real algebraic subvarieties of M", i.e., the zero sets of finite collections 

of polynomial in n real variables. 
E2.: For every linear map L : M" — ^ M, the half-plane {x G M"; L{x) > 0} belongs to S". 
Pi.: For every n > 1, the family S" is closed under boolean operations, U, n and complement. 
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P2.: If A G S*", and B G S", then Ax B G S"^+". 

P3.: If A G S*", and T : M'" ^ M" is an affine map, then r(^) G S". 

Example B.l (Semialgebraic sets). Denote by M^aig the collection of real semialgebraic sets. Thus, 
A G ^2ig if ^^'^ ^^^y if ^ is a finite union of sets, each of which is described by finitely many 
polynomial equalities and inequalities. The celebrated Tarski-Seidenberg theorem states that Saig is a 
structure. □ 

Let S be an R-structure. Then a set that belongs to one of the S"-s is called ^-definable. If A, B 
are S-definable, then a function f : A ^ B is, called ^-definable if its graph := { (a, 6) G 
Ax B; b = f{a) } is S-definable. 

Given a collection A = {An)n>i, -An C 3'(M"), we can form a new structure 8(^1), which is the 
smallest structure containing S and the sets in An- We say that S(^) is obtained from § by adjoining 
the collection A. 

Definition B.2. An M-structure is called o-minimal (order minimal) or tame if it satisfies the property 

T: Any set ^4 G is a finite union of open intervals (a, 6), —00 < a < 6 < 00, and singletons 
{r}. □ 

Example B.3. (a) (Tarski-Seidenberg) The collection M-aig of real semialgebraic sets is a tame struc- 
ture. 

(b) (A. Gabrielov, R. Hardt, H. Hironaka, (H [TTl [131) A restricted real analytic function is a func- 
tion / : M"^ — M with the property that there exists a real analytic function / defined in an open 
neighborhood U of the cube Cn '■= [—1, 1]" such that 



f{x) X eCn 

X G M" \ C„. 



we denote by Man the structure obtained from §aig by adjoining the graphs of all the restricted real 
analytic functions. Then Ran is a tame structure, and the Man-definable sets are called (globally) 
subanalytic sets. □ 

The definable sets and function of a tame structure have rather remarkable tame behavior which 
prohibits many pathologies. It is perhaps instructive to give an example of function which is not 
definable in any tame structure. For example, the function x 1— sinx is not definable in a tame 
structure because the intersection of its graph with the horizontal axis is the countable set vrZ which 
violates the tameness condition T. 

We list below some of the nice properties of the sets and function definable in a fixed tame structure 
S. Their proofs can be found in ||3]|25l. We will interchangeably refer to sets or functions definable 
in a given tame structure § as definable, constructible or tame. 

>- {Piecewise smoothness of tame functions.) Suppose A is a definable set, p is a positive integer, 
and / : A — M is a definable function. Then A can be partitioned into finitely many definable sets 
Si, . . . , Sk, such that each Si is a C^-manifold, and each of the restrictions is a C^-f unction. 
>- (Triangulability.) For every compact definable set A, and any finite collection of definable subsets 
{Si, . . . , Sk}, there exists a compact simplicial complex K, and a definable homeomorphism <I> : 
\K\ — )• A such that all the sets ^^^{Si) are unions of relative interiors of faces of K. 
^ {Dimension.) The dimension of a definable set A C M" is the supremum over all the nonnegative 
integers d such that there exists a submanifold of M" of dimension d contained in A. Then 
dim^ < 00, and dim(cl(74) \ ^) < dim^. 
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^(Definable selection.) Any tame map f : A ^ B (not necessarily continuous) admits a tame 
section, i.e., a tame map s : B ^ A such that s{b) G f~^{b), V6 G B. 

>- {Local triviality of tame maps) If f : A B is a tame continuous map, then there exists a tame 
triangulation of B such that over the relative interior of any face the map / is a locally trivial fibration. 
>- (The o-minimal Euler characteristic) There exists a function Xo : S — > Z uniquely characterized 
by the following conditions. 

. xo{x u y) = xo{x) + xo{Y) - xo{x n y), vx, y g s. 

• If X G S is compact, then Xo{X) is the usual Euler characteristic of X. 
{Finite volume.) Any compact A;-dimensional tame set has finite A;-dimensional Hausdorff measure 

{Uniform volume bounds.) If / : A ^ 5 is a proper, continuous definable map such that all the 
fibers have dimensions < k, then there exists C > such that 

'^^{r^{b)) < C, ybEB. 

a 

Suppose X is a C^, oriented Riemann manifold of dimension n. We denote by Qk{X) the space 
of /c-dimensional currents in X, i.e., the topological dual space of the space Q,^p^{X) of smooth, 
compactly supported A;-forms on X. We will denote by 

the natural pairing. The boundary of a current T G ^k{X) is the {k — l)-current defined via the 
Stokes formula 

{a,dT) ■.= {da,T), Va G 
For every a G n''{X), T G ^m{X), k < m define a n T G nm-k{X) by 

{P,anT) = {aAP,T), ypen^^+\X). 

We have 

(/3, d{a n T)) = ( dp, {a n T), ) = (a A dp, T) 
= {-l)^{d{a AP)-daA(3,T) = andT) + {-l)''+\p, da n T) 

which yields the homotopy formula 

d{ar\T) = (-i)'^*^g"(an9r- (da)nr). (b.i) 

We say that a set 5 C M" is locally subanalytic if for any p G M" we can find an open ball B 
centered at p such that S n S is globally subanalytic. 

Remark B.4. There is a rather subtle distinction between globally subanalytic and locally subanalytic 
sets. For example, the graph of the function y = sin(x) is a locally subanalytic subset of M^, but it is 
not a globally subanalytic set. Note that a compact, locally subanalytic set is globally subanalytic. □ 

If 5 C is an orientable, locally subanalytic, submanifold of of dimension k, then any 
orientation org on S determines a /c-dimensional current [S, 015] via the equality 

{a,[S,oTs]) ■■= J^a, \faen%t{R^). 

The integral in the right-hand side is well defined because any bounded, A;-dimensional globally 
subanalytic set has finite fc-dimensional Hausdorff measure. For any open, locally subanalytic subset 
J7 C we denote by [S, ots] D U the current [S D U, or^]. 
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For any locally subanalytic subset X C M" we denote by Cfc(X) the Abelian subgroup of r2fc(M") 
generated by currents of the form [S, or^], as above, where 01(5") C X. The above operation [5, or^Jn 
U, U open subanalytic extends to a morphism of Abelian groups 

efc(x) BT^Tnu € ek{x n u). 

We will refer to the elements of Cfc(X) as subanalytic (integral) k-chains in X. 
Given compact subanalytic sets ^4 C X C M" we set 

2;fc(X,yl) = {T G efc(M"); suppTcX, suppSTcA}, 

and 

'Bk{X,A) = {dT + S- TeZk+i{X,A)), SeZk{A)}. 

We set 

^kiX,A) :=Zk{X,A)/'Bk{X,A). 
R. Hardt has proved in lITTll that the assignment 

{X, A)^J{,{X,A) 

satisfies the Eilenberg-Steenrod homology axioms with Z-coefficients. This implies that !K,{X, A) 
is naturally isomorphic with the integral homology of the pair. 

To describe the intersection theory of subanalytic chains we need to recall a fundamental result of 
R. Hardt, [10, Theorem 4.3]. Suppose Eq, Ei are two oriented real Euclidean spaces of dimensions 
no and respectively ni, / : i?o — ^ is a real analytic map, and T G Qno-ciE^) a subanalytic current 
of codimension c. If y is a regular value of /, then the fiber f~^{y) is a submanifold equipped with 
a natural coorientation and thus defines a subanalytic current [f~^{y)] in Eq of codimension ni, i.e., 
^do-di i^o)- We would like to define the intersection of T and [/ ^(y)] as a subanalytic 
current (T, /, y) G S„p_c-ni (^o)- It turns out that this is possibly quite often, even in cases when y 
is not a regular value. 

Theorem B.5 (Slicing Theorem). Let Eq, Ei, T and f be as above, denote by dVE^ the Euclidean 
volume form on Ei, by a;„^ the volume of the unit ball in Ei, and set 

Jlf{T) := {y G Ei; codim{suppT)nf'^{y) > c + ni, codim(supp5T)n/"^(y) > c+m + l }. 
For every e > and y £ Eiwe define T »£ f'^{y) G ^lno~c-ni (Eq) by 

{a,T.,f-Hy)):=^j^{{rdVE,)Aa,Tn{f-\B,{y))), e Q^^r""' (Eo)- 

Then for every y G 3i/(T), the currents T »£ f~^{y) converge weakly as e > to a subanalytic 
current {T, f, y) G Sno-c-m {Eq) called the /-slice ofT over y. Moreover, the map 

JlfBy^ {T,f,y) Ged^-c-dA^") 
is continuous in the locally flat topology. □ 

Appendix C. Normal cycles of subanalytic sets 

We follow the presentation in [21]. Let V be an oriented real Euclidean vector space of dimension 
n. Denote by its dual, and by the unit sphere in V^. We identify the cotangent bundle T*V 
with the product x V. We have two canonical projections 
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Let {—,—) : V"^ X V — )• M denote the canonical pairing 

The Euclidean metric (— , — ) on V defines isometrics (the classical lowering/raising the indices op- 
erations) 

{x^,y) = {x,y), (e,y> = (e^y), yx,yeV, ^ G V^- 

Let a S $7^ (T* V) denote the canonical 1-form on the cotangent bundle. More explicitly, ifx^,...,x"' 
are Euclidean coordinates on V, and ^i, . . . , ^„ denote the induced Euclidean coordinates on V^, 
then 

a = ^idx\ 

i 

We denote by w G Q'^{T*V) the associated symplectic form 

ui = —da = dx* A d^j. 

i 

For any closed subanalytic subset X C x V we denote by Cfc(X) the Abelian group of 
subanalytic, A;-dimensional currents with support on X; see Appendix|B] If S" G Ca:(S^ x V), and 
^ G , we denote by the p-slice of S over ^, 

As explained in Appendix |Bj the slice 5^ exists for all ^ outside a codimension 1 subanalytic subset 
of and it is supported on the fiber p~^{$,) n supp S. More precisely, is well defined if the fiber 
P~^{0 ^ supp S has the expected dimension, dim S — dim S^. 

If S is the current of integration along an oriented A;-dimensional manifold, then for generic ^ the 
slice 5^ is the current of integration along the fiber S n p^^{i) equipped with a canonical orientation. 
In general, the slice gives a precise meaning as a current to the intersection of S with the fiber p~ ^ (^) , 
provided that this intersection has the "correct" dimension. 

If X C V is a compact subanalytic set, ^ G S^, and x G X we set 

Xli>i{x) ■= {y ^ X] C(y)>C(x)}, := 1 - limx(Sr(2;)nX5>g(3,)), 

where x denotes the Euler characteristic of a topological space. If x G V \ X we set ix{i, x) = 0. 
For generic ^ G S^, we have ix{i, x) = 0, for all but finitely many points x G X. 
We have the following existence and uniqueness result due to J. Fu, [7, Tm. 3.2]. 

Theorem C.l. Let X tea compact subanalytic subset ofV. Then there exists exactly one subanalytic 
current N G C„_i(S^ x V) satisfying the following conditions. 

(1) The current N is a cycle, i.e., dN = 0. 

(2) The current N has compact support. 

(3) The current N is Legendrian, i.e., 

{aUri,N) = 0, ;Vr/ G ^^""^(S'^ x V). 

(4) For any smooth function if G C°°(S^ x V) we have 
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Remark C.2, Using |l5j Thm. 4.3.2.(1)] we deduce that the equaUty (iv) is equivalent with the 
condition 

where denotes the canonical 0-dimensional current determined by the point x). The points x 
for which i{x, ^) should be viewed as critical points of the function — : X — > M; see llT4l §5.4]. 
Thus, the slice records both the collection of critical points of —^\x and their Morse indices. □ 



Definition C.3. The cycle N whose existence and uniqueness is postulated by Theorem C.l is called 



the normal cycle of the compact subanalytic set X and it is denoted by AT"^. Using the metric 
identification between the unit sphere in and the unit sphere S{V) C V we will think of as 
a (n — 1) -dimensional cycle on S{V) x V. □ 

Example C.4. (a) If X is a compact smooth submanifold of V, then can be identified with the 
integration current defined by the total space of the unit sphere bundle associated to the normal bundle 
of the embedding X ^ V. 

(b) If X is a bounded domain in V with sufficiently regular boundary dX, then we have a unit outer 
normal vector field 

n:dX^ S{V) 

and the normal cycle is the integration current defined by the graph of the above map. 

(c) If / : V — [0, c«) is a proper, C^, subanalytic function, then the normal cycle of of the sublevel 
set {/ < e} converges as e — to the normal cycle of the level set {/ = 0}. 

(d) For a very intuitive description of the normal cycle of a compact PL subset we refer to ElllTl- □ 



One can show that the map that associates to a compact subanalytic set its normal cycle is injective; 
see [JjillJ- This means that a compact subanalytic set is completely determined by its normal cycle. 
The actual reconstruction process is based on a "motivic" Radon transform. 

One very useful property of normal cycles is the inclusion-exclusion property, tli Thm. 4.2], 1211 
§4] 

for any compact subanalytic sets X, Y. 

The notion of normal cycle is closely related to the concept of curvature measure. Let us observe 
that the cotangent bundle is equipped with several canonical 5*0 (V) -invariant n-forms. To describe 
them fix an oriented orthonormal basis (ei, . . . , e„) of V . Denote by (x^, . . . , x") the associated 
Euclidean coordinates, and by ^i, . . . , x„ the dual coordinates on . For t > we define 

nv,t = {dx^ + d^i) A • • • A (dx" + td^n) G n'^iV"^ x V). 

Set 

and denote by dp the radial vector field on (V^ \ 0) x V 

j 

We set 

Sk = dpCk = dpJdsk, k = l,...,n. 
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On (F^ \ 0) X F we have 

_ n n 

^dpJ nv,t = dpJ /\ {dx^ + tdij) = dpA l\{ dx^ + td{psj) ) 
i=i i=i 

= Y,{-lf-^Sk l\ [dx^ + tdCj) = J2i-'^)''~^Sk A (dx^ + i^pdsj + tsjdp) 

k jj^k k 

We denote by r/v t the restriction of jdpJ ^v,t to x V . Along this manifold we have /? = 1 and 
we deduce 

n-l 

r]v,t = ^{-if'^Sk l\ ( dx^ + tdsj) =: ^ t^Vn~i~j 
k j^k j=0 

We denote hy u>k the volume of the unit fc-dimensional ball, and by crk-i the area of its boundary. 
Then 

k 

We set 

77fc := ??fe, A: = 0, ...,n-l 

Using the metric identification between V and we will think of the forms fjk as forms onS{V) x 
V, the unit sphere bundle associated to the tangent bundle of V. 
If X is compact subanalytic set, then the quantities 

\k{X) := (??fe,Ar^) 

are called the curvature measures of X. We have the celebrated Weyl tube formula, |[26l . 

Theorem C.5. /fX C V is a smooth compact submanifold ofV of dimension m and Tr{X) denotes 
the tube of radius r around X, 

Tr{X) := {v €V; dist(t;, X) < r } 

then 



k=0 



Moreover 

A„(X) =vol(X), Xo{X) = x{X). 

In general, Am-2fe-i(^) = 0, while Xm-2kiX) can be expressed as the integral with respect to the 
Riemannian volume on X of an universal polynomial of degree k in the curvature of the induced 
metric on X. □ 



In general, for any compact subanalytic set X we have 

Ao(X) = x{X). 



60 BRANDON ROWEKAMP 

Example C.6. Suppose dim V = 2. Using polar coordinates p, in the plane V"^ we have 

^i = /9cos0, ^2 = psin^ 

d^i = cos 9dp — p sin ddO, d^2 = sin 9dp + p cos 9d9, 
^v,t = (dx^ + t cos 6dp — tp sin 9d9) A (dx^ + t sin ^dp + tp cos 
= dx^ A dx'^ + t(cos 6*^^ - p sin 6'd6') A dx"^ + tda;^ A (sin 9dp + p cos 6'd6') + t'^pdp A d6l, 

-9pf2v,t = cos9dx^ — sm9dx^ + tpd9, 

Vv,t '■= cos 9dx^ — sm9dx^ + td9, 

rio = d9, fjo = —d9 = ^d9, 

(Tl ZTT 

o 1 1 1/ o 1\ 

r/i = cos9dx — sm9dx , fji = — r/i = -(cos^dx — sin^dx j. 

Suppose now that X is a compact smooth domain with connected boundary. We orient the boundary 
using the outer-normal-first convention and we fix an arc length parametrization of the boundary 
compatible with this orientation 

Xl=Xl{s), X2=X2{s), SG[0,L], 

where L is the length of the boundary. Denote by t{s) the unit tangent vector 

r(s) = (xi(s),X2(s)) 

and by n{s) the unit outer normal. The frame (n(s) , t(s)) and, using polar coordinates, we can write 

n{s) = (cos0(s), sin^(s) ). 
The normal cycle of X is the current of integration given by the oriented, closed path 

[0,L] 9 siA {cos9{s),s\n9{s); xi(s),X2(s)) G S{V) x V. 

Since frame {n{s),T{s)) is positively oriented we deduce that t9s) is obtained from n{s) via a 
counterclockwise rotation by | . This implies that 

x'i{s) = — sin 6'(s), X2(s) = cos 9{s). 

Hence 



{fi 



, AA^) = £ $*r/i = IJ\ {x'lisf + ) = = ^ X perimeter of X. 



□ 
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Figure 17. A 5 x 5 grid of pixels properly indexed. 



Appendix D. The Approximation Algorithm 



In this section we describe Algorithm |4.12 loosely in terms of a computer program. The input 



for this algorithm is a pixelation P£{S), where S" is a compact set. Since PeiS) is compact, it is 
contained within some m x m rectangle of e-pixels in the plane. Associate to each pixel a value 



in [1, m] n Z X [1, m] n Z indexed as a matrix (see Figure 17 1. Using this labeling we can encode 
Pe{S) as am X m matrix A where the entry aij = 1 if and only if the pixel associated to lies 
in P£{S), and is otherwise. Throughout the algorithm we will refer to the matrix A constructed in 
this manner. The notation C[i, j] will indicate the center of the pixel (in P£{S)) corresponding to the 
entry a-ij. 

The algorithm requires a choice of spread to function. Set 

a{e) = [m^'J, 
where r is a fixed rational number r G (g, 1). Note that 

lim e(cr(e))^ = oo, limecr(e) = 0. (D.l) 

e\0 e\0 

The output of the algorithm will be a PL set Ss that decomposes in a canonical fashion as a finite 
union of trapezoids with vertical bases. We will refer to such regions as poly trapezoids. We allow for 
degenerate trapezoids, such as points, segments, or triangles. 

The algorithm uses several basic subroutines. The first one is the the subroutine stack. Its input is 
a list 

C = Cl, . . . , Cm, ^1 = 0,1, 

which will be a column from A. The output of Stack is a list of nonnegative integers 

n(C); bi<ti<b2<t2<--- < bn{C) < tn(C), 

where n{C) is the number of stacks in the column encoded by C, and the location of the bottom and 
top pixel in the j-th stack is determined by the integers bj,tj. More formally 

Cfc = 1 ^=> 31 < j < n(C) : bj < k < tj. 

If C = Ci, the i-th column of A, i.e., 

Ci = Ojjl, . . . , CLi^m 
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then we will denote the output Stack(Ci) by 

fit, < < • • • < bi^m < U,ni- 
A number l<i<m — lis called a jump point if 

rii ^ rij+i. 

The next subroutine that we need is called jump. Its input is an integer k G [1, m) and the output 
is an integer jk = jijmp(A;) defined by as follows. If 

[i G [k, m) n Z; i is a jump point } = 0, 

then we set 

jump(A;) -.= 171 + 1. 

Otherwise 

jump(A;) = min{i G [k, m) n Z; i is a jump point }. 
The noise region is determined by a finite collection of intervals 

[£i,ri], . . . , [£a,ra] C [l,m] 
where the integers £k, fk are determined inductively as follows. 

Ix = max(jump(l) - 2c7(e), l), 

n = min(m,jump(l) + 2a{e)). 
Suppose that £i,ri,. . ., £j, Vj are determined. If jump(rj) > m we stop. Otherwise we set 

£j+i = max(jump(rj) - 2(7(e), l), 

rj+i = min(m, jump(rj) + 2c7(e) ) . 
The intervals [£i,ri], . . . ,[£a, ra] may not be disjoint, but their union is a disjoint union of intervals 

[ai,bi\,...,[aj,bj], 6^ < Oj+i. 
The intervals [aj, bj],l < j < J axe the noise intervals. The intervals 

[1 , ai] , [6i , 02] , . . . , [6 j-i , a j] , [bj, m] 

are the regular intervals. 

The heart of the algorithm consists of two procedures, one for dealing with the noise intervals 
and and the other for dealing with the regular intervals. These procedures will return a number of 
polytrapezoids. 

First some notation. Given a collection of points 

Bo, To,. . . , Bn, Tjv G 

such that 

x{Bi) = x{Ti), y{Bi)<y{Ti), Vz = 0,...,7V, 

x{Bj-i) <x{Bj), Vl< j<iV, 

we denote by polygon (i?o, To, . . . , B^, Tjv) the region surrounded by the simple closed PL-curve 
obtained as the union of line segments 

[-Bo, Bi],. . . , [Bn^i, Bn], 

[Bn,Tn], ■ ■ ■ , [Ti,To], [ro,i?o]- 
Note that each of the quadrilaterals Bi^iBiTITi^i is a (possibly degenerate) trapezoid with vertical 
bases. 
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Consider first the regular intervals. Given a regular interval / := [p, q] we observe that the number 

of stacks rii is independent of i G [p,q]- We denote this shared number by n = n{I). 
We construct inductively a sequence of numbers ig < ■ ■ ■ < in as follows: 

• We set io = p. 

• If q — p < 2a{e) we set iV = 1 and ii = q. 

• If io, . . . , ifc are already constructed, then, iiq — ij^ < 2a{£) v/esetN = k + 1 and ife+i = 9> 

else ifc+i = ifc + o-(e)- 
Note that if q — p > a{£), then N > 1, zq = p, in = q and 

iV = 1 if q-p<a{e). 

We have 

StaCk(Ci^) = n, bi^^^i,ti^^i, . . . ,bi^^n,tik,n- 

For j = 1, . . . , n, and = 0, . . . , we denote by B^j the center of the e-pixel corresponding 
to the element entry bi^^j in the column Ci^.. Similarly we denote by T^j the center of the pixel 
corresponding to the entry tj^, j of the column Cj^.. For 1 < j < n{I), we set 

:= polygon (5o,i, To J, . . . , BN,j,TNj). 

Define 

?(/) = U P,(/), O^reg := U P{I). 

j=l I regular interval 

Suppose now that / = [p, g] is a noise interval. We modify the column 

Cp = 1 , . . . , CLp^m 

to a column 
by setting 

^p,k '~ ' 

[o, if ELp«^,A.>o. 

We apply the subroutine Stack to the new column C^and the output is 

stack(Cp) = n(/), bi < h < ■ ■ ■ < bn < tn- 
For j = 1,. . . , n{I) we set 

Bo,j:=C\p,bj], To,j:=C\p,tj], 

Bi,j := C[q,bj], Tqj := C[q,tj], 
(recall that C[i, j] is defined as the center of the pixel associated to (i, j)). Next, for j = 1, . . . , n(/) 
we define the rectangle 

:= polygon(Eo,„ro,„5i,„rij), 

and we set 

31(1) = \J Jlj{I), noise := U R{I)- 
j=l I noise interval 

The output of the algorithm is the polytrapezoid 

•P£(^) ■= -Pregular U -Pnoise- 
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